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Preface 


For  the  reader  who  is  unfamiliar  with  viscoelasticity,  Chapter  III  presents  a  general  overview 
of  some  basic  concepts.  It  also  introduces  some  of  the  notation  used  in  the  rest  of  the  document. 
Most  of  Chapter  IV  could  also  be  considered  introductory,  since  the  presentation  doesn’t  depart 
from  well-established  theory  until  about  the  end  of  the  chapter. 

The  symbol  list  was  cumbersome  enough  that  I  included  it  as  Appendix  A,  rather  than  in  the 
prefatory  pages.  The  notation  there  appears  roughly  in  the  same  order  in  which  it  is  introduced  in 
the  text.  Most  symbols  are  also  defined  where  they  are  first  introduced  in  the  text. 

This  thesis  would  not  have  been  possible  if  1  had  not  received  plenty  of  help  and  encourage¬ 
ment.  My  faculty  advisor,  Dr.  Torvik,  showed  exceptional  patience  and  tolerance.  I  also  received 
some  timely  advise  from  Dr.  Jones  in  the  Materials  Laboratory,  which  profoundly  affected  the  final 
content  of  the  thesis.  Of  course,  my  wife  and  family  deserve  credit  for  enduring  the  long  hours, 
and  for  shoring  up  my  confidence.  My  wife  has  always  had  more  faith  in  me  than  I  have  in  myself. 

All  students  in  the  GAE90D  graduating  class  are  sobered  by  the  memory  of  a  classmate  who 
died  leist  year  in  an  automobile  accident.  Such  a  tragedy  broadens  our  perspective  on  school,  and 
on  life.  May  we  all  keep  a  keen  eye  on  our  priorities. 

In  memory  of  Wayne  P.  Wilsdon. 

Eugene  T.  Cottle 
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AFlT/GAF/i: 

I  Abstract 

'i 

tor  ninny  yrnrs,  viscoelastic  materials  have  been  used  in  damping  treatments  to  control 
vibration.  In  order  to  optimize  a  damping  treatment  for  maximum  damping,  a  design  engineer 
must  be  able  to  predict  frequencies  of  vibration  of  damped  normal  modes  and  loss  factors  as  a 
function  of  the  design  parameters.  Theoretical  models  of  layered  beams  have  been  developed,  but 
to  date,  exact  solutions  have  assumed  lioundary  conditions  to  be  the  same  for  ail  layers. 

In  this  thesis,  equatioins  of  motion  for  a  damped,  layer'd  beam  are  solved  for  the  longitudinal 
displacement  of  both  elastic  layers,  as  well  a.s  the  l.ater.al  displacement  of  the  composite  beam.  Test 
cases  we;--'  run  on  a  computer  to  verify  the  validity  of  ;Tie  equations,  and  to  comp.are  results  with 
previously  published  approximate  methods.  Results  show  that  the  software  compares  well  with 
theory  for  casc.s  in  which  simpie  theoretical  solutions  are  available. 

To  determine  the  elToct  of  mixed  boundary  conditions  on  design  optinrizaiion.  the  system 
d-amping  factor  was  calculnti’d  ,as  ;t  function  cT  ceiiiesiv.-  lay-  r  th:''kr!ess  using  Irtli  t:;o  i  rL-'.icusly 
piiblisli'.d  approximat',-  nvthod.s  and  the  i’,«  w  lV.rmu!.ati';n.  P'-r:n:tting  ir.ixe  i  1  .;:;ndar;.  '-.-n  iitiens 
tsa.s  .'liown  to  significantly  eifcct  di  .sign  (..■[itimizatiiM!. 


DAMPING  OF  LAYERED  BEAMS  WITH  MIXED  BOUNDARY  CONDITIONS 


I.  Introduction 

Engineers  frequently  face  the  problem  of  excessive  vibration  in  structures,  especially  in  the 
aerospace  industry.  One  approach  to  eliminating  or  reducing  excessive  vibration  is  to  use  materials 
which  tend  to  dampen  vibration.  Materials  can  be  designed  into  the  structure,  or  they  can  be 
added  to  an  existing  structure.  The  latter  approach  is  more  common,  since  vibration  problems 
usually  only  show  up  after  the  structure  is  put  into  service.  Damping  materials  are  commonly 
added  to  a  structure  in  the  form  of  a  damping  tape.  Damping  tape  is  a  viscoelastic  adhesive  on  a 
foil  substrate. 

A  design  engineer  is  interested  in  maximizing  system  damping  for  a  structure.  In  order  to 
maximize  damping,  it  is  necessary  to  calculate  the  damping  as  a  function  of  the  design  parameters. 
The  system  damping  is  a  function  of  the  frequency  of  vibration  of  steady  state  harmonic  motion. 
To  calculate  system  damping,  then,  it  is  necessary  to  solve  the  equations  of  motion  of  the  structure 
of  interest  for  the  damped  normal  modes  and  resonant  frequencies.  Theoretical  models  of  damped 
structures  have  been  developed  (18),  and  the  literature  suggests  that  optimum  configurations  ex¬ 
ist  (13). 

One  of  the  structures  which  has  been  investigated  extensively  is  a  layered  beam.  A  schematic 
of  a  layered  beam  is  shown  in  Figure  1.  It  consists  a  base  layer,  generally  the  structure  of  interest, 
covered  by  a  viscoeleistic  layer  and  an  elastic  constraining  layer.  The  constraining  layer  induces 
high  shear  strains  in  the  adhesive  layer.  The  shear  strain  in  the  viscoelastic  layer  is  the  primary 
damping  mechanism  (7). 

Solutions  to  the  equations  of  motion  of  the  layered  beam  have  considered  boundary  conditions 
only  on  the  composite  beam  (13).  In  other  words,  the  boundary  conditions  for  all  layers  are  the 
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Figure  1.  Constrained  Layer  Damping  Treatment 

same.  However,  in  practice,  it  is  likely  that  the  layers  of  the  damped  structure  will  not  have  the 
same  boundary  conditions.  The  question  which  needs  to  be  answered,  then,  is  how  would  the 
calculation  of  an  optimum  system  design  be  affected  by  permitting  the  boundary  conditions  of  a 
layered  beam  to  be  different  for  the  individual  layers. 

In  this  thesis,  the  equations  of  motion  are  solved  and  system  damping  are  calculated  for  a 
layered  beam  in  which  the  elastic  layers  are  permitted  to  have  different  boundary  conditions.  The 
effect  of  permitting  mixed  boundary  conditions  on  the  optimum  system  design  is  investigated. 

A  brief  historical  overview  is  first  presented  to  show  how  the  present  work  relates  to  previously 
published  theoretical  work.  Following  the  historical  overview,  the  reader  is  introduced  to  some  of  the 
basic  concepts  of  viscoelastic  damping.  The  equations  of  motion  are  then  derived  using  Hamilton’s 
principle  (QAd).  The  boundary  conditions  are  shown,  and  boundary  condition  matrices  are  shown 
for  several  important  cases. 

Numerical  results  are  presented  which  demonstrate  the  validity  of  the  equations  and  software, 
and  which  illustrate  the  effect  of  considering  mixed  boundary  conditions  on  the  optimum  design. 
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Some  of  the  factors  which  were  consiJered  in  programming  the  equations  are  discussed  briefly  Test 
cases  are  considered  for  which  simple  solutions  existed.  It  is  shown  that  the  complex  frequencies 
predicted  by  the  software  compare  well  with  the  results  predicted  by  the  simple  theory.  Complex 
frequencies  and  mode  shapes  are  shown  which  illustrate  the  effect  of  considering  mixed  boundary 
conditions.  In  addition  to  mode  shapes  of  lateral  vibration  of  the  composite  beam,  mode  shapes  of 
longitudinal  vibration  of  the  elastic  layers  are  also  shown.  Finally,  sample  optimization  curves  are 
shown.  Damping  is  calculated  as  a  function  of  adhesive  layer  thickness  for  several  cases.  It  is  shown 
that  the  formulation  presented  presented  here  predicts  an  order  of  magnitude  different  optimum 
thickness  than  do  previously  published  approximate  methods.  It  is  also  shown  that  considering 
mixed  boundary  conditions  significantly  affects  the  optimum  design. 
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II.  Historical  Perspective 


There  have  been  literally  hundreds  of  theoretical  treatments  of  viscoelastic  damping  mech¬ 
anisms  over  the  past  four  decades.  There  have  also  been  comprehensive  literature  reviews  by 
recognized  experts  in  the  field  (18).  This  perspective  will  review  only  what  might  be  considered 
milestones  in  the  theoretical  development,  and  consider  the  general  trends  in  the  research,  partic¬ 
ularly  as  it  concerns  the  subject  of  this  thesis. 

For  a  number  of  years  before  a  theoretical  treatment  was  available,  viscoelastic  material  had 
been  used  as  a  damping  treatment  in  the  form  of  damping  tape  (7).  The  damping  tape  consisted 
of  viscoelastic  material  on  a  foil  substrate. 

Governing  equations  for  structures  with  a  constrained  viscoelastic  layer  were  presented  by 
Kerwin  in  1958,  who  cited  earlier  theoretical  work  by  Oberst  and  Licnard  (7).  The  viscoelastic 
characteristics  of  a  material  were  modelled  by  permitting  the  shear  modulus  to  take  on  complex 
values.  Kerwin  identified  shear  strain  as  the  primary  damping  mechanism  in  a  layered  beam. 
DiTaranto  (5)  solved  the  sixth  order  auxiliary  equation  for  longitudinal  displacement  of  the  beam. 
Mead  and  Markus  (8)  reformulated  the  equations  in  terms  of  the  lateral  displacement  of  the  beam, 
listed  possible  boundary  conditions,  and  demonstrated  orthogonality  of  damped  normal  modes.  In 
this  early  research,  while  a  broad  category  of  boundary  conditions  were  acknowledged,  solutions 
were  restricted  to  the  simply  supported  case.  More  complicated  boundary  conditions  presented 
numerical  difficulties  (13).  It  is  also  important  to  note  that  boundary  conditions  were  considered 
on  the  composite  beam,  but  not  on  individual  elastic  layers. 

In  1978,  Trompette  et  al  (20)  presented  a  finite  element  solution  for  the  problem  of  mixed 
boundary  conditions.  At  about  the  same  time,  Rao  (13)  presented  theoretical  solutions  and  exact 
frequencies  and  loss  factors  for  all  possible  boundary  conditions  on  the  composite  beam.  Rao 
also  presented  approximate  formulae  which  satisfactorily  approximated  the  exact  solutions  for 
material  loss  factor  <  0.3. 
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During  the  past  decade,  several  important  theoretical  advan'-es  were  made.  Bagley  and 
Torvik  (2)  introduced  fractional  calculus  models  for  transient  analysis  of  viscoelastically  damped 
structures.  Simply  stated,  fractional  calculus  permits  differentiation  and  integration  of  arbitrary 
order  rather  than  just  of  integer  order.  In  other  words,  a  derivative  may  be  represented  by  d^f/dx^, 
where  A  is  not  restricted  to  integer  values.  Miles  and  Reinhall  (11)  included  effects  of  thickness 
deformation  of  the  adhesive  layer.  Finally,  Torvik  (19)  derived  the  equations  of  motion  and  bound¬ 
ary  conditions  for  a  multiple  layered  configuration  in  which  the  longitudinal  displacements  of  the 
elEistic  layers  are  considered  separately.  That  treatment  laid  the  groundwork  for  an  exact  solution 
in  which  mixed  boundary  conditions  are  permitted. 

This  thesis  uses  an  approach  similar  to  Rao  (13)  to  solve  the  equations  of  motion  presented  by 
Torvik  (19)  for  lateral  displacement  of  the  composite  beam,  as  well  as  for  longitudinal  displacements 
of  the  separate  elastic  layers.  The  approach  can  be  generalized  to  multiple  layers  as  was  suggested 
by  Torvik  (19). 
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III.  Preliminary  Concepts 


A  brief  explanation  of  viscoelasticity  will  be  presented,  to  familiarize  the  reader  with  basic 
concepts  and  notation  used  in  the  rest  of  the  document.  The  loss  factor  is  defined,  and  a  distinction 
is  made  between  the  material  loss  factor  and  the  beam  equivalent  loss  factor. 

The  stress-strain  relation  for  a  one  dimensional,  linear,  elastic  material  is  generally  represented 
by  (6;89) 

a  =  Ec  (1) 

where 

a  =  stre.$s 
€  =  strain 

E  =  Young’s  modulus,  an  experimentally  determined  propor¬ 
tionality  constant 

Material  viscoelasticity,  in  contrast  to  elasticity,  permits  the  relationship  to  include  dependence 
not  only  on  stress  and  strain,  but  on  stress  rates  and  strain  rates  (3:14-25).  For  example,  one 
simple  model  of  viscoelasticity,  called  the  Kelvin  or  Voigt  model,  for  an  isotropic  material,  could 
be  written 

cr  =  £’ie -p  £’2^  (2) 

where 

El  and  £2  =  experimentally  determined  proportionality  constants 

t  =  time 

The  dependence  of  stress  on  the  strain  rate  is  called  ‘creep’,  because,  for  a  constant  stress,  the 
strain  continues  to  change  with  time.  This  can  be  seen  by  solving  Eqn  2  for  f  eis  a  function  of 
cr  and  t.  In  a  similar  way,  strain  can  be  related  to  stress  rate.  Stress  rate  dependence  is  called 
‘relaxation’,  because  it  permits  stress  to  change  with  time  while  strain  is  held  constant. 
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A  more  general  model  could  be  expected  to  include  not  only  a  dependence  on  the  first  deriva¬ 
tives  of  stress  or  strain  with  respect  to  time,  but  on  higher  order  derivatives  as  well.  The  stress-strain 
relationship  for  a  general,  one  dimensional,  linear  viscoelastic  material  could  be  written  (3:14) 


^  dt”'  “  ^  £f<n 

m=0  n=0 


(3) 


where 

Cm  and  bn  =  emperical  proportionality  constants 

M  and  N  =  the  number  of  derivatives  and  proportionality  constants 
considered 

For  a  simple  elastic  material,  M  =  N  =  0.  In  the  Kelvin  model,  M  ~  0  and  =  1.  While  Eqn  3 
shows  dependence  on  integer  derivatives,  it  is  also  possible  to  consider  a  dependence  involving 
fractional  derivatives  (2).  In  other  words,  m  and  n  need  not  be  integers. 

For  steady  state  oscillatory  motion,  the  stress  and  strain  can  be  represented  by  periodic 
functions  of  time  (3:21): 

a  =  <Toe'"'  and  e  =  coe*"‘  (4) 


where 


(To  =  average  stress 
(0  =  average  strain 

u>  =  complex  frequency 

Substituting  equations  4  into  Eqns  3  results  in  two  complex  polynomials  in  frequency.  Factoring 
(To  and  Co  from  the  polynomials,  and  dividing  through  by  the  polynomial  on  the  left  side  of  the 
equation,  yields 


ELo^>»(»'u;) 

IIt7i=0  ‘*m(*‘*^)' 


(5) 


which  is  of  the  form 


complex  stress  =  complex  modulus  x  complex  strain 

According  to  this  argument,  it  is  possible  to  represent  the  stress-strain  relationship  of  a  general, 
one  dimensional,  linear  viscoelastic  material  with  a  single  complex  modulus,  which  can  be  experi¬ 
mentally  determined.  The  modulus  will  be  a  function  of  frequency. 

Although  Eqn  3  leads  to  an  intuitive  feel  for  meaning  of  a  complex  modulus,  it  will  be  shown 
that  most  materials  cannot  be  represented  by  a  linear  model  with  integer  derivatives.  Fortunately, 
the  concept  of  a  complex  modulus  is  still  valid  for  more  complicated  constitutive  relations.  Refer¬ 
ence  (3)  gives  a  more  rigorous  treatment. 

In  analyses  of  viscoelastic  damping  treatments,  the  shear  modulus  of  the  adhesive  layer  is 
generally  taken  to  be  complex  (7),  of  the  form  G(1  -t-  The  term  t)^  is  called  the  loss  factor  of 
the  material.  The  c.ibscript  v  will  be  used  to  differentiate  between  the  loss  factor  of  the  viscoelastic 
material,  and  the  equivalent  loss  factor  of  the  composite  beam.  The  beam  loss  factor  will  be  denoted 
Tji,  and  is  explained  below.  Note  that  is  the  ratio  of  the  imaginary  and  real  parts  of  the  complex 
modulus.  The  real  and  the  imaginary  parts  of  the  complex  modulus  are  called  the  storage  and  loss 
moduli,  respectively  (3:22). 

Although  typically  only  the  shear  modulus  of  the  adhesive  layer  is  considered  complex, 
the  elastic  modulus  may  also  be  complex.  This  concept  may  be  used  to  analyze  multi-layered 
beams  (19).  The  treatment  in  the  chapters  which  follow,  and  the  software  developed  under  this 
effort,  permit  a  complex  elastic  modulus. 
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The  loss  factors  of  important  damping  materials  have  been  experimentally  determined  and 
characterized  in  a  damping  design  guide  published  in  1985  (17).  Data  in  the  design  guide  is 
presented  in  a  graphical  format  called  a  reduced  temperature  nomogram  (17:1-13,503-511).  A 
designer  can  either  read  damping  data  from  the  nomogram,  or  calculate  it  using  an  analytical 
representation.  The  analytical  representation  is  a  curve  fit  of  the  experimental  data.  To  use  the 
analytical  representation,  the  designer  first  calculates  a  reduced  temperature,  called  qt,  then  uses 
the  reduced  temperature  to  calculate  the  complex  frequency,  from  the  following  formulae; 


logOT  = 


1 


+ 


T  A(1)J 
b 


+  2.303 


\A(1) 


-  fc  log 


T 

A(l) 


LA(1)  A(l)2 


A{4) 


{T-A{\)) 


where 


(6) 


T 

a 

b 

Ca 

Cb 

Cc 

Da 

Db 

Dc 

De 

A(1)...A(6) 


temperature 
{DbCc-CbDc)/De 
(CaDc  —  DaCc)  IDe 


[A{2)  A(1)J 

1  1 

A{2)~W) 

A(o)-A{4) 

1  11^ 

A(3)  A(l). 

1  1 

A(3)  "  A(l) 

.4(6) -A(4) 

DbCa  -  CbDa 

curve  fit  parameters  specified  for  material 
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and 


where 


G=  5(1)  + 


_ m _ 

1  +  5(5)  [5(3)/(/.)]®(«)  + 


(7) 


G 

i 


f 

5(1)... 5(6) 


Complex  shear  modulus 

frequency 

curve  fit  parameters  specified  for  material 


The  temperature  and  frequency  should  be  specified  in  units  which  are  consistent  with  the  parame¬ 
ters  yl(l) . . .  j4(6)  and  5(1) .. .  5(6).  Note  that  in  Eqn  7  the  complex  modulus  is  not  a  simple  ratio 
of  polynomials  in  frequency,  as  suggested  by  Eqn  3.  As  mentioned  earlier,  this  does  not  invalidate 
the  concept  of  a  complex  modulus.  It  only  suggests  that  most  materials  cannot  be  adequately  rep¬ 
resented  by  the  linear  constitutive  relation  in  integer  derivatives.  The  derivation  in  the  following 
chapters  permits  the  modulus  to  be  a  function  of  frequency,  but  does  not  make  any  assumptions 
regarding  the  functional  form  of  that  dependence. 

There  are  a  number  of  ways  to  quantify  the  damping  of  a  system.  They  include  the  loss 
factor,  percent  critical  damping,  viscous  damping  factor,  log  decrement,  complex  frequency,  etc. 
All  these  quantities  are  related  (15:2-7).  To  avoid  the  inevitable  confusion,  only  two  quantities  will 
be  used  in  this  thesis.  One  will  be  the  complex  damped  natural  frequency  of  vibration: 


U!  =  W/J  +  IjJli 


(8) 
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where 


ui  =  complex  frequency  of  vibration 
utji  —  real  part  of  the  complex  frequency 
ujj  =  imaginary  part  of  the  complex  frequency 

The  other  will  be  the  loss  factor  rj.  The  distinction  between  the  beam  loss  factor,  rji,  and  the 
material  loss  factor  rj^ ,  will  be  made  in  a  moment. 

The  relationship  between  the  loss  factor  and  the  complex  frequency  can  be  seen  by  solving 
the  equation  of  motion  of  a  second  order  system.  The  equation  of  motion  of  free  vibration  of  an 
undamped  harmonic  oscillator  is  (10:19) 

mu  +  ku  —  0  (9) 

where 

m  =  mass 

u  =  displacement 

k  =  spring  constant  or  elastic  modulus 

Damping  can  be  introduced  into  the  equation  by  letting  the  spring  constant  be  complex: 

mx  +  k(  \  +  iT))x  —  0  (10) 

This  is  solved  in  the  standard  fashion,  by  assuming  x  —  e"*'*,  noting  that  ui  may  be  complex.  After 
some  simple  manipulation,  it  can  be  shown  that 

(  1  +  177)  =  —  —  W;  +  i)  (11) 
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Equating  the  imaginary  parts  leads  to  a  relationship  between  the  loss  factor  and  the  complex 
frequency: 


2u>  ftujj 


(12) 


or 


2u)/ 

T]  r: -  for  uij  <C 

UlR 


(13) 


For  a  homogeneous  beam,  the  loss  factor  is  simply  the  loss  factor  of  the  material.  However, 
for  the  layered  beam,  the  loss  factor  defined  by  Eqn  12,  will,  in  general,  be  different  than  the  loss 
factor  of  the  material  in  the  adhesive  layer.  For  this  reason,  subscripts  will  be  used  to  differentiate 
between  the  two.  The  loss  factor  of  the  material  will  be  called  rjt, ,  and  the  loss  factor  of  the  beam, 
defined  by  Eqn  12,  will  be  called  rjj. 

The  preliminary  concepts  which  have  been  shown  here  are  enough  to  permit  us  now  to  derive 
the  equations  of  motion  of  a  layered  beam. 
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IV.  Derivation  of  Equations  of  Motion 


The  equations  of  motion  will  be  derived  for  a  three-layer  beam  using  Hamilton’s  princi¬ 
ple  (9:44).  To  apply  Hamilton’s  principle,  it  is  necessary  to  first  find  expressions  for  the  kinetic  and 
potential  energy.  This  derivation  is  similar  to  work  which  can  be  found  in  the  literature  (8)  (13). 
It  will  be  shown  how  the  equations  of  motion  differ  from  previous  work. 

Since  the  goal  is  to  calculate  ‘damped  normal  modes’  (8),  the  derivation  which  follows  will 
consider  only  the  homogeneous  case,  i.e.  externally  applied  loads  will  not  be  included  in  the 
potential  energy  terms. 

Figure  1  shows  the  directions  of  the  lateral  displacements,  u,  and  the  longitudinal  displace¬ 
ments,  u.  Lateral  displacements  are  perpendicular  to  the  axis  of  the  beam,  and  longitudinal  dis¬ 
placements  are  along  the  axis  of  the  beam.  Only  the  plane  strain  case  will  be  considered  (14:220). 
Plane  strain  implies  that  displacements  ‘out  of  the  page’  are  zero. 

The  kinetic  energy  of  the  beam  is 
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where 

L  =  length  of  the  beam 
6  =  width  of  the  beam 
PeO  =  density  of  base  (elcistic)  layer 
tea  =  thickness  of  the  base  layer 
HjO  =  lateral  displacement  of  the  bcise  layer 
Ujo  =  longitudinal  displacement  of  the  base  layer 
p„i  =  density  of  1st  viscoelastic  layer 
<vi  =  thickness  of  the  1st  viscoelastic  layer 
Vi,  1  =  lateral  displacement  of  the  1st  viscoelastic  layer 

Ui.i  =  longitudinal  displacement  of  the  1st  viscoelastic  layer 
Pel  =  density  of  1st  elastic  layer 
tji  =  thickness  of  the  1st  elastic  layer 
Del  =  lateral  displacement  of  the  1st  elastic  layer 
tiei  =  longitudinal  displacement  of  the  1st  elastic  layer 

Kinetic  energy  of  rotation  is  assumed  to  be  negligibly  small  in  Eqn  14.  This  is  reasonable 
when  displacements  are  small  compared  to  the  length  of  the  beam.  Eqn  14  can  be  simplified  by 
cissuming  the  and  the  kinetic  energy  due  to  the  longitudinal  displacements  is  small  compared  to 
that  due  to  lateral  displacements.  Also  assuming  v  =  Veo  =  and  eissuming  unit  width,  6=1, 
reduces  the  kinetic  energy  to 

where 

(pt)^  =  PfOteO  +  Pel<el  +  (1®) 

is  the  equivalent  areal  density  (mass  per  unit  area)  of  the  beam. 
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The  potential  energy  is 


V  = 


1 

2 


+  EeoteO 


+  Eeltei 


dx 


(17) 


where 

Eeo  =  modulus  of  elasticity  of  the  base  (elastic)  layer 
Gvi  =  shear  modulus  of  the  1st  viscoelastic  layer 
Eel  =  modulus  of  elasticity  of  the  1st  elastic  layer 

and 


7i.i 

di 

Df 


Uei 


-  UeO  ,  ,  dv 


tel  +  Go 

2t,i 

Eeoi%  Eelifi 

12(1-4)  12(1-4) 


It  has  been  assumed  (5)  (8)  that 


(18) 

(19) 

(2C) 


2 

til 


Ee\t 


t'  1 


1  /  dUei  dUeo\ 
2\dx^dx) 


'>  Ge0te07e0  +  Geltel7el 

«  +  Eeoteo{^y 


(21) 

(22) 


Eqn  21  states  that  the  shear  strain  energy  in  the  elastic  layers  is  small  compared  to  the  shear 
strain  energy  of  the  adhesive  layer.  This  is  recisonable  since  the  damping  treatment  is  inducing 
large  shear  strains  in  the  adhesive  layer.  Eqn  22  states  that  the  strain  energy  in  the  adhesive  layer 
due  to  longitudinal  displacement  is  small  compared  to  the  longitudinal  strain  energies  of  the  elastic 
layers.  This  is  reasonable  since  the  longitudinal  strain  of  the  adhesive  layer  cannot  be  greater  than 
the  strains  in  the  elastic  layers,  and  the  elastic  modulus  of  the  adhesive  layer  will  be  small  compared 
to  the  el^lstic  layers. 
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Hamilton’s  principle  (9:44)  requires: 


where  S  is  the  variational  operator  ( 14:442-443) . 

The  equations  of  motion  and  boundary  conditions  are  obtained  by  integrating  Eqn  23  by 
parts.  Tlie  integration  by  parts  is  carried  out  in  such  a  way  that  terms  involving  derivatives  of 
Si:,  SUf()  and  Sud  are  eliminated  from  the  integral,  and  only  terms  involving  Sv,  Sii(0  and  Su^i 
themselves  remain  under  the  integral.  In  doing  so,  terms  are  extracted  from  the  integral  which 
contain  conditions  which  must  be  satisfied  at  ti  and  t2,  and  at  x  =  0  and  L.  The  variations  are  =  0 
at  and  t2.  This  is  equivalent  to  saying  that  initial  conditions  will  be  specified.  The  conditions 
which  must  be  satisfied  at  x  =  0  and  L  are  the  boundary  conditions.  They  will  be  listed  later. 
What  results  from  the  application  of  Hamilton's  principle  and  the  integration  by  parts,  then,  is 
an  integral  expre.ssion  which  must  be  satisfied,  and  a  set  of  boundary  conditions.  The  integral 
expression  which  must  be  satisfied  is 


Sv-Dt—,S,: 


'^EeoteO  'n-OUeo  +  Ee\te\  "" 

dx^  dr- 


(24) 


The  integral  must  be  identically  0  for  all  x.  The  variations  6v,  6uco  and  are  permitted  to 
vary  arbitrarily  and  inuependently  over  the  domain  (14:445)  from  x  =  0  to  i  =  L,  so  the  coefficients 
of  the  variations  inside  the  integral  must  all  vanish.  This  leads  to  the  equations  of  motion: 


^  d^v  ,  ,  5^1)  ,  r  1  ( duei  dUfo\  ,  d~v 

{Pt)r  (^—  _  _  j  +  d, _ 

^  f  Ucl  -  UtO  ,  j  dv 

+  dr  ^ 

r.  ,  ^  f  «el  -  «f0  ,  J  dv 

V  "“i  ^ 


The  boundary  conditions  were  extracted  from  Eqn  23  during  the  integration  by  parts,  as 
conditions  which  must  be  satisfied  at  x  =  0  or  x  =  L.  They  are 

UeO  =  0 

Uel  =  0 

(26) 

v  =  0 

dv/dx  =  0 

The  first  of  Eqns  26  states  that,  for  the  base  layer,  either  the  longitudinal  displacement  must 
be  zero,  or  the  ‘slope’  of  the  longitudinal  displacement  must  be  zero.  Requiring  the  slope  to  be  zero 
is  equivalent  to  requiring  the  end  to  be  force-free.  Longitudinal  displacement  must  be  unrestrained. 
The  second  equation  requires  the  same  to  be  true  for  the  constraining  layer.  The  third  “quation 
requires  either  the  lateral  displacement  of  all  layers  to  be  zero  or  be  free  of  shear.  It  must  be  either 
restrained  laterally,  or  free  to  displace  laterally.  Note  that  the  expression  for  the  shear  force  at  the 
end  is  somewhat  more  complicated  than  comparable  expressions  for  homogeneous  beams  (9:134). 
It  could  be  considered  an  equivalent  shear  force  for  the  composite  beam.  The  last  equation  requires 


- 

0 

or 

- 

0 

or 

-  +  Ui 

1 

bv 

dx 

= 

0 

or 

d-v/dx~  —  0  or 
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either  the  slope  to  be  zero  or  the  moment  to  be  zero.  It  must  be  either  prevented  from  rotating, 
or  free  to  rote*'-  These  equations  must  be  satisfied  at  both  ends  of  the  beam,  so  they  actually 
represent  a  set  of  eight  boundary  conditions. 

These  equations  of  motion  and  boundary  conditions  compare  with  published  results  in  the 
early  literature  (5)  (8)  (13),  except  Eqns  26  consider  the  longitudinal  displacements  of  the  elastic 
layers  separately  (19).  Note  that  in  Eqns  25,  Uco  and  u^x  appear  only  as  a  difference.  Considering 
boundary  conditions  on  the  composite  beam  is  equivalent  to  considering  boundary  conditions  on 
the  difference  between  and  Ud.  In  that  case,  only  six  equations  result  (13). 
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V.  General  Solution  of  the  Equations  of  Motion 


In  the  literature  (8)  (13)  (19),  Eqns  25  are  usually  combined  into  a  single  sixth  order  equation. 
In  order  to  consider  mixed  boundary  conditions,  it  will  be  necessary  to  obtain  expressions  explicitly 
for  Ueo  and  Uji.  This  is  more  easily  done  if  the  three  separate  equations  are  solved  simultaneously. 
It  is  possible  to  obtain  the  same  results  by  solving  the  single  sixth  order  equation  for  v,  and 
substituting  the  resulting  expression  back  into  Eqns  25  to  solve  for  u^o  and  Uei,  but  that  procedure 
is  more  tedious. 

Since  the  problem  of  interest  is  steady  state  oscillatory  motion,  displacements  will  be  assumed 
to  be  periodic  functions  of  time.  It  is  customary  (7)  (5)  (8)  to  assume  a  solution  of  the  form 

V  -  e'‘^‘Aef^ 

Ueo  =  (27) 

Uei  = 

Note  that  this  ^lssumes  the  solution  is  separable.  In  other  words,  it  assumes  that  the  displacements 
can  be  written  as  a  function  of  x  alone  times  a  function  of  t  alone.  Assuming  a  separable  solution 
implies  that  mode  shapes  are  expected  which  are  not  a  function  of  time,  as  in  the  steady  state. 
The  orthogonality  of  the  solutions  (8)  guarantees  the  uniqueness  of  the  steady  state  solution,  but 
a  transient  solution  of  this  form  may  not  be  unique. 
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Substituting  Eqns  27  into  Eqns  25  and  dividing  through  by  e‘“'  and  results  in  a  set  of 
simultaneous  algebraic  equations.  They  can  be  written  in  matrix  form  as  follows 


[DtP^  -  Gt,idip 

Gvldip  [£^eO<eOP^  —  (Gvi/<„i)] 


—Gvid\p 

Gv\/t^\ 


—G^\d\p 


Gv\/ix}\ 


\Ec\te\p'^  — 


'  •< 

'  \ 

A 

0 

X  < 

B 

.  =  , 

0 

C 

0 

(28) 


Finding  the  solution  to  Eqns  28  is  an  eigenvalue  problem,  although  it  is  not  posed  in  the  same 
general  form  as  most  eigenvalue  problems  (9:138).  Still,  it  can  be  solved  in  the  same  way.  A 
nontrivial  solution  requires  the  determinant  of  the  three  by  three  matrix  =  0; 


0 


_ (pOt 

Dt 


<.1^1  ,  1 

[  Dt  <t,i 


1 


EelU 


,  {pt)ru,-G.i  /  I 

^  Drtvl  \Ee}te\ 


Ef  teO 

+ 


EeoteO 


(29) 


This  is  the  characteristic  equation.  It  differs  from  the  characteristic  equation  published  in  the 
literature  (5)  (13)  in  that  it  has  a  double  zero  root.  The  double  zero  root  is  the  result  of  considering 
longitudinal  displacements  of  the  individual  elastic  layers  separately.  Note  that  only  even  powers 
of  p  appear  in  the  equation.  This  means  that  the  roots  will  come  in  pairs.  For  p  ^  0,  the  roots  are 
founu  by  first  solving  for  the  roots  of  a  cubic  equation  in  p^,  then  taking  the  positive  and  negative 
square  root  of  each.  The  result  will  be  six  complex  values  In  positive  and  negative  pairs. 

For  each  value  of  p  which  satisfies  Eqn  29,  Eqns  28  are  linearly  dependent.  For  each  p  91^  0, 
the  solution  is  determined  to  within  one  arbitrary  constant.  For  p  ~  0,  the  solution  is  determined 
to  within  two  arbitrary  constants. 
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The  pj  are  eigenvalues  of  the  system.  Mead  and  Markus  (8)  demonstrated  that,  for  pj  ^  0, 
the  associated  eigenfunctions  are  orthogonal.  For  p  —  0,  the  eigenfunctions  can  be  made  orthogonal 
by  choosing  proper  linear  combinations  of  them  (9:142).  The  expansion  theorem  then  guarantees 
that  any  function  satisfying  the  homogeneous  boundary  conditions  can  be  represented  as  a  linear 
combination  of  the  eigenfunctions  (9:143).  This  means  that,  for  solutions  of  the  form  of  Eqns  27, 
the  mode  shapes  will  be  uniquely  determined  by  a  linear  combination  of  the  eigenfunctions. 

For  p  ^  0,  the  eigenfunctions  are  determined  by  solving  for  B  and  C  in  terms  of  A  in  Eqns  28. 
Since  the  equations  are  now  linearly  dependent,  any  two  of  the  equations  may  be  used.  Using  the 
first  two  of  Ens  28  results  in  the  following  relation: 


and 


{pt)j.u^-p]DT 


c,  = 


^  P^Dt  -p]G^iUxd\  - 


d\Pj 


Aj  =y^Aj 


If  the  second  and  third  of  Eqns  28  are  used,  the  relations  become 


(30) 


(31) 


B,= 


—  EeltelGyltyldlPj 


A,^I3,A, 


(32) 


and 


C, 


EeOtcoGy\tyld\Pj 


\_BeQitoEe\tc\ty\Pj  —  G„i  {EeoteO  +  E'eitel)  j 


7^.4; 


(33) 


Eqns  30  and  31  are  equivalent  to  Eqns  32  and  33  for  pj  which  satisfy  Eqn  29.  These  expressions 
relate  the  lateral  displacement  of  the  composite  beam  to  the  longitudinal  displacements  of  the 
individual  elastic  layers. 


21 


For  p  ^  0,  there  are  now  only  six  unspecified  constants,  A\  ...As-  Since  there  are  eight 
boundary  conditions,  two  more  constants  are  expected.  The  two  additional  constants  will  result 
from  the  double  zero  root  in  the  characteristic  equation. 

For  p  =  0,  the  solution  is  of  the  form: 

i;  =  (At  +  Asx) 

ti,o  =  (Bt  +  Bsx) 

=  e‘“‘(C7+Csx)  (34) 


Substituting  these  into  Eqns  28  makes  the  equations  linearly  dependent. 


The  two  independant 


equations  which  remain  are 


{-^7  +  A^x)  —  Cjt,\diCg  +  G,,id\Bs  —  0 

^(Ct  +  Csx- St- BsJ-)  +  6V,d,yl8  =  0 

HI 


[35) 


Collecting  the  coefficients  of  x'^  and  x'  and  setting  them  =  0  yields 


—  At  —  GxiidiCs  +  GfidiB^  —  0 


- —  (C't  —  Bt)  +  GvidiA^  —  0 

Hi 

-(Pt)u)^As  =  0 


Gx,\ 

ivl 


(Q  -  Sg) 


0 


(3G) 


For  w  0 


Aj  =  As  =  0 

St  =  Ct 
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B6  =  Cs 

For  u!  =  0,  A7  and  .dg  are  unspecified  (corresponding  to  rigid  body  lateral  translation  and  rotation), 
and 


Bs  =  Cg 

Bt  —  €7  =  dit^iAs 


If  u.'  =  0  satisfies  the  boundary  conditions,  that  case  must  be  considered  separately.  It  will  not  be 
treated  here. 


The  general  solution  to  Eqns  25  is  then 


V  =  € 


UeO  =  e 


U,1  =  e 


iwt 


J  =  l 
6 

{pj  aj€^’^) + B7  +  R^x 

6 

+  B7  +  Bsx 

j  =  i 


(37) 


The  solution  is  now  determined  to  within  eight  arbitrary  constants,  Ai  S7  and  Bs- 

This  was  to  be  expected,  since  there  are  eight  boundary  conditions  which  need  to  be  satisfied. 
There  is  also  one  unspecified  parameter,  the  complex  frequency  of  vibration,  w.  It  now  remains 
to  find  the  frequency  which  satisfies  a  particular  set  of  boundary  conditions.  Specific  boundary 
condition  matrices  are  considered  next. 
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VI.  Boundary  Condition  Matrices 


Substituting  Eqns  37  into  an  appropriate  combination  of  Eqns  26  results  in  a  set  of  eight 
simultaneous  equations  in  the  eight  constants  Ai  ...ylg,  B-;  and  and  frequency.  In  matrix 
form,  it  can  be  represented  as 

[boundary  condition  matrix]  x  {column  vector  of  coefficients)  =  {zero  vector) 

Nontrivial  solutions  will  exist  only  when  the  determinant  of  the  boundary  condition  matrix  is 
zero.  Solving  the  equations  involves  finding  the  complex  frequency  which  drives  the  determinant  to 
zero.  In  contrast  to  the  boundary  conditions  on  the  composite  beam,  where  the  boundary  condition 
matrix  is  six  by  six  (13),  for  mixed  boundary  conditions,  in  general  the  matrix  will  be  eight  by  eight. 
However,  it  will  be  shown  that  the  matrix  reduces  to  a  six  by  six  when,  at  both  ends  of  the  beam, 
the  boundary  conditions  are  the  same  for  each  elastic  layer.  In  this  chapter,  boundary  condition 
matrices  are  presented  for  several  important  Ccuses.  Boundary  conditions  were  selected  which  would 
be  most  illustrative  of  the  effects  of  considering  the  elastic  layers  separately.  A  summary  of  the 
boundary  conditions  considered  is  shown  in  Table  1.  The.^e  boundary  conditions  are  used  in  the 
numerical  results  presented  in  Chapter  VI 1. 


Both  Base  and  Constraining  Layers  Fixed  at  x  =  0,  Free  ai  x  =  L 


The  case  where  both  the  base  and  constraining  layers  fixed  at  x  =  0  and  free  at  x  =  L  is 
of  interest  because  it  duplicates  solutions  wind,  are  available  in  the  literature  (16:5-24).  It  will  be 


Table  1.  Summary  of  Boundary  Conditions  Considered 


Base  Layer 

Restraining  Layer 

X  0 

X  ~  L 

X  =  0 

X  —  L 

fixed 

free 

fixed 

free 

fixed 

free 

iree 

free 

pinned 

pinned 

free 

free 

pinned 

pinned 

rigidly  connected 

free 

used  for  comparison  purposes.  The  boundary  conditions  are;  at  x  =  0 


dv/dx  =  0  zero  slope 

u  =  0  zero  lateral  displacement 

Ueo  =  0  zero  longitudinal  displacement  of  base  layer 

Uei  =  0  zero  longitudinal  displacement  of  constraining  layer 


and,  at  X  =  L 


duca/dx  =  0  zero  longitudinal  force  on  base  layer 

duei/dx  =  0  zero  longitudinal  force  on  constraining  layer 

d'v/dx~  =  0  zero  moment  on  composite  beam 

^  j  ,  [uei  -  u,o  ,  j  ^ 

=  0  zero  shear  on  composite  beam 
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Substituting  Eqns  37  into  Eqns  38  and  39  produces  a  set  of  eight  simultaneous  equations 
in  the  eight  constants,  Ai, . . .  ,Ae,  B-  and  Bs-  In  this  case,  the  constants  By  and  Bg  can  be 
eliminated  from  the  system,  reducing  the  number  of  simultaneous  equations  to  be  satisfied  to  six. 
The  equations  are,  in  matrix  form 


1 

1 

1 

1 

1 

1 

Pi 

P2 

P3 

P4 

P5 

P6 

^2 

'I'a 

’J'4 

^5 

'i'6 

^2P2«''‘“ 

^4P4eP*^ 

'J'ePee^*^ 

p|eP=^ 

p2eP=^ 

ple^D 

ple^<-^ 

/ 

’ 

Ai 

0 

A2 

0 

Az 

y  ~  4 

0 

A4 

0 

A5 

0 

Ae 

\  J 

0 

where 


(40) 


'^}  =  P}  -  T> 


and 


dihj  -  p,)  A  d\i^iPj 


Dtp] 


Base  Layer  Fixed  at  x  -  0  and  Free  at  x  =  L;  Constraining  Layer  Free  at  x  =  0  and  at  x  =  L 

This  is  an  example  of  mixed  boundary  conditions.  Comparing  this  case  to  the  previous  one 
will  illustrate  the  effect  of  considering  mixed  boundary  conditions  on  the  loss  factor.  The  boundary 


26 


conditions  are;  at  x  =  0 


dvjdx  =  0  zero  slope 

t;  =  0  zero  lateral  displacement 

Ueo  =  0  r.ero  longitudinal  displacement  of  base  layer 

dxie\ldx  —  0  zero  longitudinal  force  on  constraining  layer 


and,  at  X  =  L 


dueo/dx  =  0  zero  longitudinal  force  on  base  layer 
dufi/dx  =  0  zero  longitudinal  force  on  constraining  layer 
d'^v/dx‘  =  0  iero  moment  on  composite  beam 

=  0  zero  shear  on  composite  beam 
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Since  the  boundary  conditions  on  the  elastic  layers  are  different  at  i  =  0,  the  constants  By 
and  Bg  cannot  be  eliminated.  The  equations  are,  in  matrix  form 

1  1  1  1  1  10  0 
Pi  P2  P3  P4  P5  P6  0  0 

01  02  03  04  03  /?6  1  0 

7lPl  72P2  73P3  74P4  75P5  76P6  0  1 

0iPieP'^  02P2e^'‘^  03P3e^^^  04P4e^*^  06Pe.e^^^  0  1 

7iPieP‘^  72P2e^’^  73P3e'’^^  75P5e’’*^  76)P6e'’‘'^  0  1 

p\e^'^  p\ef^^  p\e^^^  p\e^*^  p\e^^^  0  0 

hie^'^  h.2e^^^  0  0 


Although  this  cannot  be  reduced  to  six  simultaneous  equations,  the  sparsity  of  the  last  two  rows 
simplifies  the  calculation  of  the  determinant.  It  is  apparent  that,  at  the  most,  three  six  by  six 
determinants  need  to  be  calculated. 
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Base  Layer  Pinned  ai  x  =  0  and  at  x  =  L;  Constraining  Layer  Free  at  x  =  0  and  Free  at  x  =  L 

This  boundary  condition  was  selected  for  comparison  with  known  solutions.  It  is  closest  to  the 
simply  supported  composite  boundary  condition  which  was  the  first  case  solved  in  the  literature  (8). 
The  boundary  conditions  are;  at  both  1  =  0  and  x  =  L 


V 

=  0 

zero 

d^v/dx^ 

=  0 

zero 

=  0 

zero 

dufi/dx 

=  0 

zero 

lateral  displacement 
moment  on  the  composite  beam 
longitudinal  displacement  of  bcise  layer 
longitudinal  force  on  constraining  layer 


(44) 
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The  equations  are 


1 

1 

1 

1 

1 

1 

0 

0 

0 

PI 

P2 

vl 

vl 

vl 

vl 

0 

0 

/?! 

/?2 

03 

A 

03 

06 

1 

0 

7lPl 

72P2 

73P3 

74P4 

75P5 

76P6 

0 

1 

03ef^^ 

0^eP*^ 

066’’^^ 

1 

L 

7iPieP‘^ 

72P2e^^^ 

73P3eP'’^ 

74P4eP*^ 

75P5e’’*^ 

76P6eP*^ 

0 

1 

p\eP^^ 

pleP^^ 

p\eP*^ 

p^eP^i 

ple”-^ 

0 

0 

fPit. 

gPsL 

gPtL 

gP^L 

gPiL 

0 

0 

•V 

Ai 

0 

A  2 

0 

A3 

0 

A4 

►  =  t 

0 

0 

els 

0 

B- 

0 

Ss 

0 

(45) 


In  this  case,  by  subtracting  row  five  from  row  three,  87  can  be  eliminated  and  the  number  of 
simultaneous  equations  reduced  to  seven. 


Base  Layer  Pinned  at  x  ~  0  and  at  x  =  L;  Constraining  Layer  Rotates  with  Base  Layer  at  x  =  0 
aiid  Free  at  x  =  L 

This  is  a  variation  on  the  simply  supported  case  in  the  previous  section,  and  will  be  used 
to  illustrate  the  effect  of  permitting  mixed  boundary  conditions.  The  constrainig  layer  is  ‘rigidly 
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connected'  to  the  base  layer  at  x  =  0,  such  that  all  layers  rotate  together  at  that  end.  The 
constraining  layer  will  have  a  non-zero  enforced  displacement. 

At  first  glance  this  appears  to  be  a  non-homogeneous  boundary  condition.  All  the  boundary 
conditions  considered  so  far  have  been  homogeneous,  i.e.  all  zeros  appear  on  the  right  hand  side 
of  the  equations.  An  arbitrary  enforced  nonzero  displacement  is  a  non-homogeneous  boundary 
condition.  However,  the  enforced  displacement  in  this  case  is  not  arbitrary,  but  is  related  to  the 
slope  of  the  beam.  For  small  rotations,  at  x  =  0 


/  tco 

V  2 


dv 

dx 


ttj\d 


dv 


(46) 


Since  dv/di  contains  only  terms  in  Ai  . . .  Ae,  this  displacement  can  be  included  in  the  left-hand 
side  of  the  boundary  condition  equation,  and  the  boundary  condition  is  seen  to  be  homogeneous. 
The  boundary  conditions  are;  at  x  =  0 


dh-/dx^- 

titO 

Uei 


0  zero  lateral  displacement 

0  zero  moment  on  the  composite  beam 

0  zero  longitudinal  displacement  of  base  layer 

(dv/dx)  constraining  layer  rotates  with  the  beise  layer 


and  X  =  h 


(47) 


V  = 

d'v/dx'^  = 

«e0  = 

duti/dx  = 


0  zero  lateral  displacement 
0  zero  moment  on  the  composite  beam 
0  zero  longitudinal  displacement  of  base  layer 
0  zero  longitudinal  force  on  constraining  layer 


(48) 
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The  matrix  equation  is 


1 

1 

1 

1 

1 

1 

0 

0 

Pi 

P2 

0 

P3 

P4 

Ps 

Pe 

0 

0 

/?1 

02 

03 

A 

05 

06 

1 

0 

Fi 

F2 

Fa 

F. 

Fs 

Fe 

1 

0 

/?3eP>^ 

AeP*^ 

/?5eP»^ 

06e’’^^ 

1 

L 

7iPieP‘^ 

72P2eP=^ 

73P3e’’^^ 

74P4eP*^ 

75P5e^‘^ 

76P6eP^^ 

0 

1 

p2eP3^ 

Pae'”^ 

p\eP*^ 

0 

0 

gPii 

gP^i. 

gP3t 

gPiL 

ePti- 

ep(.i^ 

0 

0 

' 

0 

>42 

0 

>43 

0 

>44 

►  —  ^ 

0 

>45 

0 

>46 

0 

Bt 

0 

Bs 

0 

>  4 

where 

Fj  =  Tj  +  U\diPj 


(49) 


Neither  nor  Bs  can  be  eliminated  in  this  case. 

In  Chapter  VII,  numerical  results  are  given  for  the  boundary  conditions  which  have  been 
shown  here. 


32 


VII.  Numerical  Results 


The  equations  presented  in  th’  previous  chapters  were  programmed  and  results  calculated 
for  some  example  cases.  The  name  of  the  program  used  in  the  calculations  is  ‘BEAM’.  The  cases 
considered  were  intended  to  do  two  things:  validate  the  equations  and  software,  and  illustrate  the 
effects  of  considering  mixed  boundary  conditions. 

The  first  set  of  cases  verify  that  the  equations  are  .^alid  and  the  software  is  working  correctly. 
This  is  done  by  demonstrating  that  the  solutions  degenerate  into  known  solutions  in  special  cases. 
The  validity  of  the  equations  and  software  are  seen  in  the  numerical  results,  as  well  as  in  the  mode 
shapes  predicted  by  BEAM.  It  will  be  shown  that  the  software  performs  as  expected. 

The  second  set  of  test  cases  demonstrate  the  effect  of  permitting  mixed  boundary  conditions, 
and  show  how  the  computed  results  compare  with  previously  published  approximate  methods. 

For  input  data,  all  cases  use  some  variation  of  the  physical  parameters  in  an  example  calcu¬ 
lation  in  a  design  guide  published  in  1985  (16:5-26).  This  was  done  to  permit  comparison  with  the 
results  in  the  design  guide.  The  data  is  shown  in  Table  2.  The  physical  properties  of  the  viscoelastic 
material  come  from  the  third  volume  of  the  design  guide,  reference  (17:263-270). 

Computer  Programming  Approach 

Computer  software  was  written  to  implement  the  equations  shown  in  the  previous  sections. 
The  main  program  is  called  BEAM.  In  this  section  the  programming  approach  is  discussed.  The 
numerical  techniques  are  explained  briefly.  A  detailed  description  of  input  data  is  included  as 
Appendix  B.  Programming  notes  are  in  Appendix  C.  Appendix  C  also  describes  the  computer  and 
compiler  used,  data  availability,  etc. 

As  pointed  out  by  Rao  (13),  solving  the  equations  for  boundary  conditions  other  than  simply 
supported  presents  some  numerical  difficulty.  The  boundary  condition  matrix  tends  to  be  ill- 
conditioned  (4:175).  This  can  be  attributed  primarily  to  the  exponential  terms.  The  roots  of 
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Table  2.  Basic  Data  Used  as  Starting  Point  for  All  Test  C2ises  (16:5-26)  (17:263-270) 

1 .  Base  Layer,  aluminum 

2.  Adhesive  Layer,  SoundCoat  DIAD  no.  606 

3.  Constraining  Layer,  aluminum 


Variable 

Value 

Units 

L 

0.254 

meter 

EeO 

6.8900el0 

Pascal 

teO 

5.0800e-3 

meter 

PeO 

2.7700e  3 

(kg)/meter^ 

t'eO 

0.3300 

<1,1 

2.5400e-4 

meter 

Pvl 

9.6900e  2 

(kg)/meter^ 

Eel 

6.8900el0 

Pascal 

<el 

2.5400e-4 

meter 

Pel 

1.0 

(kg)/meter^ 

I'el 

0.3300 

A(l) 

335.0 

degrees  K 

A(2) 

280.0 

degrees  K 

A(3) 

390.0 

degrees  K 

A(4) 

0.7e-l 

degrees  K/Hz 

A(5) 

0.1142 

degrees  K/Hz 

A(6) 

0.3e-l 

degrees  K/Hz 

B(l) 

0.2e6 

Pciscal 

B(2) 

1200.0e6 

Pascal 

B(3) 

0.3e7 

Hz 

B(4) 

0.55 

B(5) 

1.5 

B(6) 

0.1 
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the  characteristic  equation  come  in  pairs.  Even  modest  values  of  pj  produce  both  very  large  and 
very  small  terms  in  the  matrix.  In  BEAM,  two  things  are  done  to  ensure  accurate  calculation 
of  the  determinant.  First,  double  precision  arithmetic  is  used  in  all  calculations.  Second,  total 
pivoting  is  used  in  the  gaussian  elimination  algorithm,  rather  than  scaled  partial  pivoting  (4:159). 
The  determinant  is  calculated  using  triangular  factorization  (4:160).  Triangular  factorization  is 
a  gaussian  elimination  algorithm.  Rows  and/or  columns  are  exchanged  during  the  elimination  to 
reduce  round-off  error.  In  scaled  partial  pivoting,  only  rows  are  interchanged.  In  total  pivoting, 
both  rows  and  columns  are  interchanged  as  necessary.  Total  pivoting  is  slower  than  scaled  partial 
pivoting  because  it  involves  more  overhead.  It  is  also  more  complicated  because  the  program  must 
keep  track  of  both  row  and  column  interchanges.  For  these  reasons,  total  pivoting  is  seldom  used. 
In  this  case,  total  pivoting  was  necessary  to  ensure  accuracy. 

While  it  would  clearly  be  possible  to  locate  zeros  of  the  determinant  using  a  root  finding 
algorithm  such  as  Muller’s  metljod  (4;]20),  there  are  disadvantages  to  using  this  approach.  Most 
root  finders  (4:72-127)  assume  prior  knowledge  of  approximate  locations  of  zeros,  and  they  may 
not  converge  if  they  are  not  in  the  vicinity  of  a  zero.  Also,  they  cannot  be  expected  to  find  all 
zeros.  The  magnitude  of  the  boundary  condition  determinant  tends  to  increase  as  w  increases,  so 
if  good  initial  estimates  arc  not  available  a  root  finder  is  not  likely  to  converge  to  a  root. 

In  BE.4M.  zeros  were  located  using  a  simple  grid  search  rather  than  a  more  complex  root  find¬ 
ing  algorithm.  The  grid  search  is  described  below.  A  more  efficient  routine  should  be  implemented 
In  practice,  but  for  the  purposes  of  investigating  the  locations  of  the  minima,  the  grid  search  weis 
very  useful,  and  it  avoided  potential  convergence  problems.  The  grid  search  is  not  fast,  but  it  is 
robust  and  dependable.  It  is  more  practical  than  it  is  elegant. 

In  the  grid  search,  the  user  specifies  a  range  of  values  for  real  and  imaginary  frequency,  and 
an  increment  for  both  real  and  imaginary  values.  The  software  calculates  values  for  the  magnitude 
of  the  determinant  at  every  value  of  real  and  imaginary  frequency  specified.  Thi.s  is  the  grid.  The 
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software  then  checks  all  points  on  the  grid  for  relative  minima.  If  a  minimum  is  found,  the  grid 
is  refined  in  the  vicinity  of  the  minimum  to  improve  the  accuracy.  The  surface  defined  by  the 
magnitude  of  the  boundary  condition  determinant  is  relatively  smooth,  making  it  easy  to  locate 
minima.  If  the  initial  grid  is  fine  enough,  the  grid  search  will  find  all  minima  in  a  given  range. 

After  the  minima  are  found,  it  then  remains  to  determine  which  are  actually  solutions  which 
satisfy  the  boundary  conditions,  i.e.  which  are  ‘zeros’.  There  are  several  ways  to  validate  a  zero. 
First,  the  value  is  substituted  back  into  the  boundary  condition  equations.  If  the  error  in  a  boundary 
condition  equation  is  roughly  the  same  order  of  magnitude  of  the  individual  terms  of  the  equation, 
then  that  boundary  condition  is  not  satisfied  (4:169-176).  If  the  error  in  a  boundary  condition  is 
small  compared  to  the  individual  terms,  then  the  equation  is  satisfied.  Zeros  can  also  be  verified 
by  solving  for  Aj,  p-,  Bj  and  Bg,  0j  and  7j  at  that  value  of  frequency,  then  using  these  terms 
to  calculate  a  mode  shape.  The  mode  shapes  show  clearly  whether  the  boundary  condition  has 
been  satisfied.  Zeros  can  also  be  identified  by  comparison  with  previously  found  zeros.  The  actual 
solutions  tend  to  lie  on  a  smooth  curve.  Any  minima  which  lie  well  outside  the  trend  are  not  likely 
to  satisfy  the  boundary  conditions. 

Test  Cases  for  Validating  Software 

Three  test  ceises  were  considered  for  which  comparisons  could  be  made  with  theoretical  results. 
The  first  two  cases  were  for  a  cantilever  beam  with  the  constraining  layer  fixed  at  x  =  0.  This  is  the 
first  boundary  condition  listed  in  Table  1.  These  two  cases  used  input  data  for  the  layered  beam 
for  which  the  frequencies  could  be  expected  to  approach  the  frequencies  of  an  undamped  cantilever 
beam.  The  third  uses  the  third  boundary  condition  in  Table  1,  and  the  input  data  in  Table  2.  This 
boundary  condition  should  closely  resemble  the  simply  supported  case  in  the  literature. 

Case  1  -  Bare  Beam  The  first  test  ra.se  approximated  a  bare  cantilever  beam,  by  letting 
the  shear  modulus  of  the  damping  layer  be  small,  and  the  elastic  modulus  and  density  of  the 
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Table  3.  Data  Used  to  Calculate  Theoretical  Frequencies  for  Cases  1  and  2 


Variable 

Case  1 

Case  2 

Units 

E 

6.8900el0 

6.8900el0 

Pciscal 

1 

1.0925e-8 

1.2990e-8 

meter'' 

L 

0.254 

0.254 

meter 

m/ L 

14.0716 

15.0213 

(kg)/meter 

constraining  layer  be  small.  The  problem  formulation  does  not  permit  the  numbers  to  be  set  to 
exactly  zero.  In  some  cases,  setting  a  value  to  zero  would  result  in  division  by  zero.  In  other 
cases  it  would  cause  excessive  roundoff  errors.  Trial  and  error  was  used  to  find  the  smallest  values 
the  software  would  accept.  The  real  part  of  the  resulting  complex  frequencies  which  satisfy  the 
boundary  conditions  should  approach  the  natural  frequencies  of  a  simple  cantilever  beam.  The 
imaginary  part  of  the  complex  frequencies  should  approach  zero,  corresponding  to  zero  damping. 

The  theoretical  natural  frequencies  are  roots  to  the  transcendental  equation  (10:224-227) 


cos \L  cosh \L  =  — 1 


(50) 


where 

A"*  =  — 

El 

Table  3  lists  the  data  used  to  calculate  the  theoretical  frequencies,  and  Table  4  lists  the  input 
data  for  calculating  the  natural  frequencies  in  BE.4M.  Data  not  listed  in  Table  4  was  unchanged 
from  that  shown  in  Table  2.  Note  that,  in  this  case,  the  shear  modulus,  G^i,  was  held  constant  for 
all  frequencies,  rather  than  letting  it  vary  according  to  Eqns  6  and  7.  Table  5  lists  the  resulting 
natural  frequencies  for  the  first  six  modes  as  predicted  by  Eqn  50,  and  the  complex  frequencies 
predicted  by  BEAM. 

Note  in  Table  5  that  the  imaginary  frequencies  predicted  by  BEAM  are  zero  or  very  small, 
suggesting  small  damping,  as  expected.  All  real  frequencies  are  5.9%  high.  This  could  be  the  result 
of  letting  the  parameters  of  the  top  layers  be  small,  but  not  zero. 
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Table  4.  Input  Data  for  Program  BEAM  for  Cases  1  and  2 


Variable 

Bare  Beam 
Caise  1 

Rigid  Spacer 
Case  2 

Units 

Gvi 

1.0000  -h  l.Oj 

S.OOOOe  9  +  l.Oi 

Pciscal 

Pvl 

1.0000 

1.0000 

(kg)/meter^ 

Eel 

l.OOOOe  4  +  O.Oi 

6.8900el0  +  O.Oi 

Pascal 

2.5400e-4 

2.5400e-4 

meter 

Pel 

1.0 

2.7700e  3 

(kg)/meteP 

Table  5.  Case  1  -  Comparison  of  Bare  Beam  Theoretical  Frequencies  with  Values  Predicted  by 
BEAM 


Frequencies  are  in  radians/second 


Mode 

Expected 

BEAM 

Real  Imag 

error  in 
Real  freq 

1 

398.59 

422.23 

0.0000 

5.9  % 

2 

2497.92 

2646.10 

0.0000 

5.9  % 

3 

6994.25 

7409.17 

0.0000 

5.9  % 

4 

13705.93 

14519.03 

0.0000 

5,9  % 

5 

22656.89 

24000.99 

0.0001 

5.9  % 

6 

33854.47 

35853.33 

0.0001 

5.9  % 

Case  2  -  Rigid  Spacer  Case  2  approximates  an  undamped  beam  with  a  rigid  spacer.  In 
BEAM,  the  real  part  of  the  shear  modulus  of  the  viscoelastic  layer  was  made  large,  and  the  imagi¬ 
nary  part  small.  As  in  case  1,  the  theoretical  natural  frequencies  of  the  beam  were  found  by  solving 
Eqn  50.  A  beam  with  a  rigid  spacer  will  have  a  slightly  larger  moment  of  inertia  than  a  bare  beam, 
due  to  the  separation  between  the  two  eleistic  layers.  The  mass  per  unit  length  will  be  larger  than 
the  bare  beam  because  of  the  addition  of  the  mass  of  the  constraining  layer.  The  data  used  to 
calculate  the  theoretical  natural  frequencies  for  this  case  is  shown  in  Table  3.  The  input  data  used 
in  program  BEAM  for  case  2  is  shown  in  Table  4.  As  in  the  first  Ccise,  the  complex  shear  modulus  is 
not  permitted  to  vary  with  frequency  according  to  Eqns  6  and  7,  but  remains  fixed.  The  theoretical 
natural  frequencies,  and  the  complex  frequencies  predicted  by  BEAM,  are  listed  in  Table  6. 

Again,  the  imaginary  components  of  the  frequency  are  zero,  as  expected.  The  real  components 
are  all  slightly  high. 
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Table  6.  Case  2  -  Comparison  of  Rigid  Spacer  Theoretical  Frequencies  with  Values  Predicted  by 
BEAM 


Mode 

Expected 

BEAM 

Real  Imag 

error  in 
Real  freq 

1 

420.67 

445.42 

0.0000 

5.9% 

2 

2636.31 

2791.29 

0.0000 

5.9% 

3 

7381.75 

7815.16 

0.0000 

5.9% 

4 

14465.29 

15313.04 

0.0000 

5.9% 

5 

23912.16 

25310.22 

0.0000 

5.9% 

6 

35720.63 

37803.00 

0.0000 

5.9% 

It  should  be  noted  that  the  derivation  of  Eqn  50  assumes  that  shear  strain  energy  is  negli¬ 
gible  (10:221).  Permitting  the  layered  beam  to  have  a  rigid  spacer,  i.e.  very  high  shear  modulus, 
could  intuitively  be  expected  to  result  in  higher  natural  frequencies  than  the  bare  beam. 


Case  3  -  Simply  Supported  Historically,  the  first  boundary  condition  solved  for  the  layered 


beam  was  the  simply  supported  case  (8).  The  solution  is  of  the  form 


t)  =  £>sin 


(51) 


where 

D  =  a  constant 

n  =  the  mode  number 

TT  =  3.14159... 

X  =  distance  along  the  length  of  the  beam 
L  =  the  length  of  the  beam 

Although  the  literature  does  not  address  the  longitudinal  displacement  of  the  elastic  layers, 
intuitively  one  would  expect  the  displacement  of  the  constraining  layer  to  have  the  form 


Ue\  =  F  cos 


(52) 
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where 


f  =  a  constant 

This  can  justified  with  a  simple  argument.  When  the  composite  beam  is  deflected  upward  in  the 
center,  the  longitudinal  displacement  of  the  top  layer  should  be  zero  in  the  center.  It  should  also 
be  some  negative  value  at  x  =  0  and  the  same  positive  value  at  x  =  L. 

It  is  not  immediately  apparent  what  form  the  longitudinal  deflection  of  the  base  layer  would 
take,  since  it  is  restrained  at  both  ends. 

In  this  section,  it  will  be  demonstrated  that  the  BEAM  results  satisfy  Eqns  51  and  52  when 
the  bcise  layer  is  simply  supported  and  the  top  layer  is  completely  unrestrained. 

Taking  the  real  part  of  the  expressions  for  v  and  ttej  in  Eqns  25  results  in,  for  a  given  p, 

Re(i')  =  [Afl  cos(p/x)  —  A/  sin(p/x)]  (53) 

and 

Re(uei)  =  [(hfiAn  -  y/A/)cos(p/x)  -  (t^A/  +  7/A/j) sin(p;x)]  (54) 

The  subscript  R  denotes  the  real  part  and  the  subscript  I  denotes  the  imaginary  part  of  a  quantity. 
Several  conditions  must  be  satisfied  for  Eqns  53  and  54  to  reduce  to  the  form  of  Eqns  51  and  52. 
The  pj  come  in  pairs,  and  they  are  eigenvalues  of  the  system,  so  there  should  be  exactly  one  pair 
of  Pj  such  that 

PR  =  0  and  Pi  =  (55) 
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If  the  pair  of  roots  which  meet  the  condition  of  Eqn  55  are  arbitrarily  designated  pi  and  p2,  it  can 
be  shown  that  the  coefficients  j4i  . .  and  71  . .  -  76  must  satisfy  the  following  conditions: 

Ai,r  +  A2,r  =  0 
Aij  —  A2,i  ^  0 

^3,R  +  A4^r  =  0 

(56) 

>13,/  -  A41  =  0 

As.r  +  A6,R  =  0 
A5,i  -  A^  I  =  0 

ilhRAl.R  -  lljAij)  +  (-t2.RA2,R  - 'r2.lA2,l)  ^  0 

ill, rAi, I  -  71,1  Ai^r)  -  (y2.RA2.I  -  17,1  A2,r)  =  0 

(73,R-43,/J  -  73,/^3,/)  +  (74,fi^4,fi  -  74,/>1i,/)  =  0 

(57) 

(73,R>l3,/  -  73,/^3,fi)  -  (74,fi^4./  -  74,/^4,fl)  =  0 

(75,fi^5,fi  -  75,/^5,/)  +  (76,rt-46,fl  -  76,/>f6,/)  =  0 

(75,fi-45,;  -  75,;^5,h)  -  il6, rAq,]  -  16,/A6,r)  =  0 

If  the  conditions  in  Eqns  56  are  satisfied,  only  a  sine  term  will  remain  in  the  expression  for  v. 
Similarly,  if  the  conditions  in  Eqns  57  are  satisfied,  only  a  cosine  term  will  remain  in  the  expression 
for  U(i . 

Note  that  the  subscripts  in  Eqns  55  . . .  57  are  arbitrary.  The  only  requirement  is  that  the 
nonzero  terms  in  Eqns  56  and  57  must  be  associated  with  the  two  pj  which  satisfy  Eqns  55.  Any  two, 
but  only  two,  pj  are  expected  to  satisfy  Eqns  55.  In  the  numerical  results,  simply  eis  a  consequence 
of  the  way  the  software  was  written,  these  roots  will  have  ‘adjacent’  indices.  In  other  words,  they 
could  be  satisfied  by  pi  and  p2,  by  pj  and  p4,  or  by  ps  and  pe-  In  Eqns  56  and  57  the  nonzero 
terms  would  have  the  same  subscripts  eis  the  pair  of  pj  which  satisfy  Eqns  55. 
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Table  7.  BEAM  Output,  Test  Case  3,  Mode  1 


At  Omega  =  1221.86  +14.46i  Radians  per  second; 

P  Beta  Gamma 


J  Real  Imag 

1  21.15  8.49 

2  -21.15  -8.49 

3  0.00  -12.38 

4  0.00  12.38 

5  12.18  .23 

6  -12.18  -.23 

Constcuits  A.  .  . 


Real 

-.1667E-01 

0.1667E-01 

0.2589E-03 

-.2589E-03 

0.2075E-02 

-.2075E-02 


Imag 

-.8685E-02 

0.8685E-02 

-.1332E-02 

0.1332E-02 

-.5778E-03 

0.5778E-03 


Real  Imag 

0.3334E+00  0.1737E+00 
-.3334E+00  -.1737E+00 
-.5178E-02  0.2664E-01 
0.5178E-02  -.2664E-01 
-.4149E-01  0.1156E-01 
0.4149E-01  -.1156E-01 


T  Real 

1  0.69418126E-05 

2  -.82489694E-03 

3  0.25344642E+00 

4  -.25552087E+00 

5  0. 12232727E-03 

6  0.27700760E-02 

7  -.18594580E-02 

8  0. 14641401E-01 


Imag 

O.OOOOOOOOE+00 

0.12445740E-02 

0.65964374E+00 

-.65962757E+00 

-.61542895E-04 

-.11991999E-02 

0.31897464E-03 

-.25116113E-02 


It  now  remains  to  demonstrate  that  the  numerical  results  from  BEAM  satisfy  the  conditions 
specified.  In  the  test  case,  the  input  data  listed  in  Table  2  was  used.  The  shear  modulus  of  the 
adhesive  layer  was  permitted  to  vary  according  to  Eqns  6  and  7. 

For  the  first  mode,  using  the  input  data  in  Table  2,  the  requirement  on  p/  is 


TT 


VI  = 


TT 

0.254 


12.37 


The  output  of  BEAM  for  the  first  mode  is  shown  in  Table  7.  Table  8  shows  the  values  of  the 
quantities  defined  by  Eqns  56  and  57.  In  the  program  output  shown  in  Table  7,  the  constants  S7 
and  Bg  are  shown  as  the  seventh  and  eighth  values  in  the  vector  A.  This  was  done  for  convenience 
in  programming. 

In  Table  7,  pa  and  p4  satisfy  Eqns  55.  Of  the  terms  defined  by  Eqn  56,  the  ‘nonzero’  term 
is  three  orders  of  magnitude  larger  than  the  ‘zero’  terms.  Of  the  terms  defined  by  Eqn  57,  the 
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Table  8.  Conditions  on  Constants  A  and  -y  -  Case  3  Mode  1 


+  A2,/i 
Ai,/  —  A^j 
Az,r  +  Ai^R 
As, I  —  Aij 
Af, 

,R  +  Ae.R 
As,!  —  As,! 


-8.1796£’-04 

-1.2446jE;-03 

-2.0744£'-03 

+1.3190 

+2.8924f:  -  03 

+1.1377£:-03 


-  7l,/>ll,/)  +  (72,fl^2,R  -  72,/>l2,/)  = 

(7l,fi^l,/  -  7l,/>ll,ii)  -  (72,fi>l2,/  -  72,/^2,fl)  = 

(73,R>l3,fi  -  ~1a,!Az,l)  +  {~i^,RAi,  R  -  y4,!A^j)  = 

(73.fl-^3,  f  -  'ra,lA3,R)  -  {'y4,RAi,!  -  74,/>'l,fi)  = 

{ys,RAs,R  -  IS.lAsj)  +  (76,fl-A6,fi  -  76,7^6,/)  = 

(75,fi^5,/  -  1^,1  As,r)  -  {ys,RAs,I  -  76,/>l6,fi)  = 


+4.9347T'-  04 
+5.5696£'-  04 
-3.7779f:-  02 
+5.5178£'-  05 
+9.6719T'-  05 
+  1.8887E-  05 


‘nonzero’  term  is  two  orders  of  magnitude  larger  than  the  ‘zero’  terms.  While  this  difference  in 
magnitude  is  technically  large  enough  to  meet  the  requirement,  it  is  not  clearly  convincing. 


The  data  for  the  second  mode  is  much  better.  It  is  shown  in  Tables  9  and  10.  For  the  second 
mode  the  requirement  on  p/  is 


27r 

0.254 


24.74 

! 


In  Table  9  it  can  be  seen  that  ps  and  ps  satisfy  the  specified  conditions.  In  Table  10,  the  constants 
associated  with  ps  and  ps  satisfy  Eqns  56  and  57. 

The  difference  in  the  quality  of  the  calculated  results  between  the  first  and  second  modes  may 
be  due  to  the  accuracy  with  which  the  frequencies  were  located.  Both  sets  of  data  were  calculated 
in  the  same  run,  using  the  same  tolerance.  The  software  uses  several  tolerances  to  decide  how 
closely  to  refine  the  accuracy  of  a  frequency.  The  magnitude  of  the  determinant,  or  the  change  in 
magnitude  of  the  determinant  from  one  iteration  to  the  next,  or  the  difference  in  frequencies  from 
one  iteration  to  the  next,  are  all  criteria  which  could  terminate  the  search.  These  are  explained 
in  more  detail  in  Appendix  B.  One  of  these  criteria  may  have  been  satisfied  with  a  less  accurate 
calculation  of  the  first  frequency  than  the  second.  Still,  for  both  modes,  it  is  safe  to  say  that  the 
BEAM  results  degenerate  into  the  theoretical  solutions  shown  in  Eqns  51  and  52. 
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Table  9.  BEAM  Output,  Test  Case  3,  Mode  2 


At  Omega  =  4820.67  +75.75i  Radians  per  second: 

P  Beta  Gcunma 


J  Real  Imag 

1  29.95  12.00 

2  -29.95  -12.00 

3  24.19  1.00 

4  -24.19  -1.00 

5  0.00  -24.74 

6  0.00  24.74 

Constcints  A .  .  . 


Real 

-.2025E-01 

0.2025E-01 

0.3957E-02 

-.3957E-02 

0.7089E-03 

-.7089E-03 


Imag 

-. 1904E-01 
0.1904E-01 
-.3047E-02 
0.3047E-02 
-.2110E-02 
0.2110E-02 


Real  Imag 

0.4049E+00  0.3808E+00 
-.4049E+00  -.3808E+00 
-.7913E-01  0.6094E-01 
0.7913E-01  -.6094E-01 
-.1418E-01  0.4221E-01 
0.1418E-01  -.4221E-01 


J 

1 

2 

3 

4 

5 

6 

7 

8 


Real 

0.80797525E-16 
-.86216012E-11 
0.13165025E-14 
-. 11847728E-09 
0.5723148SE+00 
-.57231485E+00 
0.94137268E-03 
0.17361670E-11 


Imag 

O.OOOOOOOOE+00 
0.21418468E-10 
0.54727301E-15 
-. 18512119E-09 
-.41527191E+00 
0.41527191E+00 
0.30044278E-02 
0.25122977E-11 


Table  10.  Conditions  on  Constants  A  and  7  -  Case  3  Mode  2 


^1,7  -  7I2,/  = 

^3.fi  +  = 

^3,7  -  A4J  = 

+  Ae,;?  = 

As, 7  —  Ac,/  = 


-8.C215E'-  12 
-2.1418E;-  11 
-1,1848£’-  10 
+  1.8512E'-  10 
-1.3428E'-  10 
-0.8305 


(71,7?Ai,/?  -  71,7Ai,7)  +  ('12,rA2,R  -72.1^2,1) 
(71,7?Ai,/  -  7i,7Ai,p)  -  (72,rA;,7  -  72,1  A2^r) 
(73,/?A3,R  -  73.1  Azj)  +  (74,rA4 ,7?  -  j) 

(73,7?A3,7  -  73, 7 A3, 7?)  “  (74,7?A4.7  -  74,7 A4  /{) 
(75.RA5  R  -  75,1  A^j)  +  (76, 77 Ae, 77  -  7^,1  A^j) 
(75,7?A5,/  -  75,1  A5,r)  -  (7e,RAe,i  —  76,1  A^  n) 


+  1.1646E'-  11 
+  1.1955£;-  11 
-2.0656f;-  11 
+2.1869E’-  11 
+  1.8827E-  02 
+8,0772£'-  12 
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_  Mode  1 

_  Mode  2 

_  Mode  3 

-2.00 

Figure  2,  Nornializcd  Lateral  Displacement  -  Base  Layer  Pinnetl-Pimied,  Constraining  Layer 
Free- Free 

The  validity  of  the  BEAM  results  can  also  he  ilernonstrated  by  looking  at  the  predicttnl  mode 
shape's.  TliC  mod':'  sha;-es  <ai  calculated  Ijy  .sul)stituting  the  v.alues  of  ,  Aj.  ,  -/j.  B-  aiui 
/J,  into  L'.ins  27  (Without  the  t'*'  t'-nu).  \'alue.s  of  di.sphaceiueui  are  calculated  ;ia<i  pl.-.tted  for 
reprc.-eutativ  vali.  .s  of  r.  Note  that,  of  the  eight  constants,  one  is  still  uns|iecirn  d.  1  his  can 
Iji'  seh'ctrd  art.'itrarily,  .:uid  is  usu.ally  used  to  tiormalize  the  mode  shapes,  1  he  uK.ale  sh.ap'  s  still 
accurat.'ly  reflect  the  relative  magnitud'.'S  of  f.  t/..()  and  Idguri.'s  2,  3,  and  ■!  show  the  first 

three  mo'h’  shapes  for  tli"  dat.i  presentc'd  in  this  section.  '1  In'  luoile  shapc-s  are  normalized  with 
re;  p>''  - 1  I .  "r  .:  -  ..i''' -me!;: .  1  !i.-i,t  j.,.  t !i.-  m.a.ximuu;  lal i-r.a  1  •  h.-. pla-'-  i ! c  ut  in  Fi'eiU'"  'J 

!'  -e'  -  d  ir  pi  ace;;;' 'Ht  s  a.  and  are  adjusted  a'-'s -rdingK'.  1  he  mode  sloii'es  af'' 


eousisteiit  wi’h  tin'  s'ated  boundarv  conditions. 


Figure  3.  Normalized  Longitudinal  Displacement  of  Base  Layer  -  Base  Layer  Pinned-Pinned,  Con¬ 
straining  Layer  Free-Free 

Test  Cases  Illustrating  the  Effect  of  Considering  Mixed  Boundary  Conditions 

In  this  section,  boundary  conditions  listed  in  Table  1  are  contrasted,  to  illustrate  the  effects 
of  considering  mixed  boundary  conditions.  Numerical  results  are  also  compared  to  predictions  of 
FRED'S',  a  program  publislied  in  a  damping  design  guide  (16:5-29).  PREDY  uses  approximate 
formulai;  to  estimate  the  complex  frecpiencit's  of  a  three  layer  cantilever  beam. 


Comparison  of  Complex  Frequencies  for  Cantileier  Boundary  Conditions  This  case  compares 
the  first  and  second  boundary  conditions  listed  in  Table  1.  which  will  be  called  the  cantilever 
I  '-undary  r(jndit:ons.  '.vith  results  predieteei  :  y  PRED5  i  l6:5-'_’Ui.  The  complex  frequencie-.s  and 
loss  factors  predicted  by  both  BE.\.M  and  PRED'\'  are  shown,  as  are  mode  shapes  predicted  by 


BE.A.M. 


.0500 


-  Moae  1 

-  Mode  2 

-  Mode  3 


-.0500  -■ 

Figure  ■}.  Normalized  Longitudinal  Displacement  of  Constraining  Layer  -  Base  Layer  Pinned- 
Pmned,  Constraining  I.ayer  Free-Free 

Complex  frequencies  were  calculated  using  the  input  data  in  Table  2  exactly  as  shown.  BEAM 
calculates  real  and  imaginary  frequency  components,  and  PRED^'  calculates  real  freqtiency  and 
loss  factor.  The  qii.antitifs  are  related,  as  shown  in  Eqn  12.  Intuitively,  one  would  expect  the 
■  -ping  to  be  liigheT  when  tlie  constraining  layer  is  fixed,  because  shear  strains  are  higher. 

The  complex  frequeiicic.s  predicted  by  both  BE.AM  and  PRED'i’  are  shown  graphically  in 
Figure  o.  Table  11  .'^hows  the  numerical  data  used  in  Figure  5.  The  complex  frequencies  for  the 
two  ca.'ses  prcflicted  by  BF.N.M  are  quite  clo.se,  while  the  PRED’i’  calculations  tend  to  predict  lower 
imaginary  fr'-qui  nc].  ■,  in  if.'  I.m;:-  r  mC’d-  s.  .\il  thre*'  '-alctdat ions  :-;;ow  dirf-  rem  !■•■:im  !•  -s  f  -  -:  r-. 
It  is  inter' stmg  that  ali  tlir'  ’’  'Si-'-  s  sb'ow  tin'  m.aximum  loss  facto.r  at  a  di.Ter''!U  ni'.ei'’.  PPLD’t’ 
[iredicts  maximum  dampung  in  the  second  mode,  while  BFAM  predicts  maximum  damping  for  the 
fourth  and  fifth  mod>  s  wh'-ii  the  cottstr.ainiiig  l.ayt  r  is  fixed-free  and  free-free,  respectively.  Except 
for  the  first  mode,  fixing  tin'  constraitiing  layer  at  =  0  re.siilts  in  higher  datnpung.  as  expectc'd. 


Table  11.  Values  of  Complex  Frequency  Predicted  by  BEAM  and  PREDY  for  Cantilever  Bound¬ 


ary  Conditions 


Mode 

BEAM 

PREDY 
Real  Imag 

Vb 

Const.  Layer  Fixed- Free 
Real  Imag  rjj 

Const.  Layer  Free- Free 
Real  Imag  r/t 

1 

437.14 

2.78 

0.0127 

421.02 

4.17 

0.0198 

410.48 

4.17 

0.0203 

2 

2697.53 

36.32 

0.0269 

2678.24 

25.19 

0.0188 

2511.78 

40.46 

0.0322 

3 

7495.15 

114.71 

0.0306 

7474.49 

96.74 

0.0259 

6946.62 

105.29 

0.0303 

4 

14588.14 

234.25 

0.0321 

14570.69 

212.39 

0.0292 

13525.28 

180.31 

0.0267 

5 

23990.36 

381.06 

0.0318 

23976.29 

358.02 

0.0299 

22271.59 

260.66 

0.0234 

6 

35692.67 

547.31 

0.0307 

35681.46 

524.36 

0.0294 

33184.20 

344.23 

0.0207 

Table  12.  Complex  Frequencies  Predicted  by  BEAM  for  Simply  Supported  Boundary  Conditions 


Mode 

Const.  Layer  Rigid-Free 
Real  Imag  rn, 

Const.  Layer  Free-Free 
Real  Imag 

1 

1199.22 

14.07 

0.0235  1 

1221.86 

14.46 

0.0237 

2 

4780.24 

66.14 

0.0277 

4820.67 

75.75 

0.0314 

3 

10704.32 

161.38 

0.0302 

10753.82 

182.75 

0.0340 

4 

18942.83 

288.94 

0.0305 

18996.25 

319.32 

0.0336 

5 

29488.11 

440.83 

0.0299 

29542.05 

478.50 

0.0324 

6 

42331.97 

608.42 

0.0287 

42385.70 

651.11 

0.0307 

The  loss  factors  for  the  first  mode  are  counter-intuitive,  but  not  unreasonable.  The  damping  is  a 
result  primarily  of  shear  stra’’"  in  the  adhesive  layer,  and  the  shear  strain  is  related  to  the  difference 
in  longitudinal  displacement  of  the  elastic  layers,  «eo  —  Wei-  When  both  elastic  layers  are  fixed  at 
1  =  0,  the  shear  strain  at  that  end  is  zero. 

The  mode  shapes  predicted  by  BEAM  for  the  two  cantilever  boundary  conditions  for  this 
example  are  shown  in  Figures  6  through  11. 


Comparison  of  Complex  Frequencies  for  Simply  Supported  Boundary  Conditions  In  this  sec¬ 
tion,  the  complex  frequencies  are  shown  for  the  two  simply  supported  boundary  conditions,  which 
are  the  third  and  fourth  in  Table  1.  Again,  intuitively,  one  would  expect  the  loss  factor  to  be 
higher  for  the  case  where  the  constraining  layer  is  fixed  at  one  end,  due  to  higher  shear  strains. 
The  locations  of  the  complex  frequencies  are  shown  in  Figure  12.  The  data  used  in  Figure  12  is 
listed  in  Table  12. 
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Imaginary  Frequency,  rad/sec 


ue0 


Figure  7.  Normalized  Longitudinal  Displacement  of  Base  Layer,  All  Layers  Fixed-Free 


F'igiirc  9.  Normalized  Lateral  DisplacenieiU,  Base  Layer  Fixed  Free,  Constraining  Layer  Free  Free 


The  maximum  dam[)im;  occurs  when  the  constraining  layer  is  free  at  both  emis,  in  the  third 
mode.  Frir  ail  modi  s.  '.lamping  is  higher  when  the  constraining  layer  is  free  at  both  ends,  than 
wiii.n  it  is  constrained  to  rotat.'  with  tlie  base  layer  at  x  =  0.  Again,  although  this  riin.s  counter 
to  the  imtiai  expectations,  it  cati  be  explained  using  the  same  arguments  used  in  tlie  previous 
'"t I'.'ii .  Datt'.r'iii'g  is  'iri'.en  1  y  .nhear  strain  in  the  adhesive  layer,  which  is  related  to  the  i.hnerenie 
in  lungit udinal  ilisplaei  itient  of  the-  ehislic  layers,  n,,o  — “ci-  Hequiriiig  the  top  elastic  hiyer  lij  rotate 
witli  the  l  ase  la>er  elf-  ctivoly  iiti[<oses  a  zero  shear  strain  lioiindary  condition. 

11, '■  r  I’a"  l.r.-t  .'imply  ■-iipi'.orted  cas--  are  sho-wn  in  Figures  2  ilu'  ii-g!:  !.  '!  1  ■■ 

:i.  .  'ii.ap- -  r  ih--  •  ii  i  'impl>  •'Upport.  ■!  ca.'e  an- sliown  in  Figures  12  throiigli  la 

Ejjift  ‘if  I’t  nnittvxq  Mtrtii  Boundary  Conditions  on  Optimum  Di  styn  As  pointed  out  earlie-r, 
the  fgo.'il  c-f  the  de.sigti  engineer  is  to  tnaximize  the  loss  factor,  for  optitnum  sy.'lein  datiiping.  'lln- 
e'.f-'-t  I  if  nr,  :■;>  d  F-  aueiary  --i.in'lit  i-  -rt.-i  on  optimum  design  is  of  primary  coneern  In  this  se-'t  loii.  .-iti 


Figure  10.  Normalized  Longitudinal  Displacement  of  Base  Layer,  Base  Layer  Fixed  Free,  Con¬ 
straining  Layer  Free  Free 


exampl'’  of  an  optimum  design  calculation  is  shown.  Only  one  variable  is  considered,  the  advesive 
layer  thickness.  The  loss  factor  is  plotted  as  a  function  of  adhesive  layer  thickness  for  five  cp'^es. 
I  he  first  three  cases  consider  loss  factors  for  the  cantilever  boundary  conditions,  as  predicted  by 
FRFD\'  and  BEA.M.  They  are  shown  in  Figures  10  through  IS.  The  fourth  and  fifth  cases  are 
K.'sults  pr-'dicted  by  BLA.M  tor  the  two  siiiijdy  .■supported  boundary  conditions.  They  arc  shown  in 
f  igures  19  and  20.  The  (luantity  ot  interest,  is  the  thickness  at  which  the  loss  factor  is  a  maximum. 
1  hese  values  are  shown  in  Tabh’s  13  and  14. 


fe.r  the  ■eult  ih'V.T  t  .  illi  i.irv  e 
t  ilickne.'-ses.  .-MT-r  dei-re::.-iii_e  •i.e 
factor  would  continue  to  incre.ast'. 
of  magnitude  .as  the  length  i u'  tie- 


i  r.dit,  1:^.  ill  il.e  hiL'li'T  I-.',  ih'-  I'urvi'S  sliow  a  peak  for  small 
less  ;  ict,.,r  tiicii  incrc-.-u-cs,  .-gain  for  barge  t  liicknesses.  The  lo.s.s 
but  ftir  magnitudes  of  thickness  which  are  on  the  same  order 
le’am  .\t  these  values,  the  ritary  inertia  terms  are  no  longer 


negligii.le,  so  the  equations  cjf  motion  ar'-  not  valid  !T0;221i 


2.00  -I 
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Figure  13.  Normalized  Lateral  Displacement,  Base  Layer  Simply  Supported,  Constraining  Layer 
n.igid-Frce 


For  the  cantilever  boundary  conditions,  all  predictions  show  that  there  is  no  optimum  for 
the  first  mode,  and  BEAM  predicts  no  optimum  for  the  first  three  modes.  For  the  fourth  through 
the  sixth  modes.  BEAM  predicts  the  optimum  at  an  order  of  magnitude  larger  thickness  than 
does  PRED^^  For  the  BEAM  predictions,  if  the  constraining  layer  is  fixt;rl  at  r  =  0,  tlie  optimum 
thickn'’.ss  is  .smaller  than  it  would  be  if  the  constraining  layer  is  fixed.  Figure  21  shows  curves  for 
both  boundary  conditions  for  the  fourth  mode,  to  show  the  general  trend.  This  shows  that  fi.xing 
the  con.straining  layer  at  j  =  0  not  only  results  in  a  peak  at  ,a  smaller  adli  .'sive  layer  thickness,  but 
t  li”  pe.-ik  I"  .'d'')  higher. 

Ft  il.e  simply  sii[;pf'rted  cases,  again  ihere  i.s  no  ch-.ar  optimum  :'  t  'i,*'  !  \ver  mod'S.  In  the 
modes  where  there  is  a  peak,  for  some  modes  the  different  boundn'-  -  'ondition  h;ts  an  effect,  in 


other  modes  it  does  not. 
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Figure  M.  Normalized  Longitudinal  Displacement  of  Base  Layer,  Base  Layer  Simply  Supported, 
Constraining  Layer  Rigid-Free 


Table  13.  Approximate  Advesive  Layer  Thicknesses  for  Optimum  Damping  -  Cantilever  Boundary 
Conditions  (meters) 


Mode 

PREDV 

BEAM 

Diif.  between 

const,  layer  fixed  free  const,  layer  free  free 

cols  2  and  3 

1 

no  maximum 

no  maximum 

no  maximum 

2 

0.00020S9 

no  maximum 

no  maximum 

3 

0.000 109G 

no  maximum 

no  maximum 

4 

0.0000759 

0.0O03S0 

0.000480 

26.3  % 

5 

0.0000575 

0.000255 

0.000290 

13.7  % 

'3 

1 _ 

0.0000479 

0.000105 

0.000220 

12. S  % 

Tiijh'  I  t,  .\pproxim.ate  Adve.sive  Layer  Tliickne.^.qes  for  Optimum  Damping  -  SirnvTv  Sunp'U'tfti 
H'Uinilary  Condition.'^  (meters) 


.Mode 

eonst.  layer  free  free 

const,  layer  ristid  free 

Ditf. 

1  • 

no  maximum 

no  maximum 

2 

no  maximum 

no  maximum 

3 

0.00045 

0.00040 

11  7c 

4 

0.0002S 

0.0002S 

0  Sf, 

5 

0.00025 

0.00020 

20  7c 

6 

O.OOOIS 

0.00016 

11  7c 

\\\ 

\\\\ 


Figure  15.  Normalized  Longitudinal  Displacement  of  Constraining  Layer,  Base  Layer  Simply  Sup 
ported,  Constraining  Layer  Rigid-Free 
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Fizure  It'i.  Hf'am  L'j.'S  Factor  vs.  Thirkness  as  Predicted  hv  .Ml  Lavers  Fi.xi'd-Frc’e 


System  Damping  Factor 


Tfiickness  -  meters 


Figure  IV.  Beam  Loss  Factor  vs.  Thickness  as  Predicted  by  BEAM  -  All  Layers  Fixed-Frei 
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Fidurn  1-^.  Ream  Loss  Factor  vs.  Thickness  a-s  Predictcfd  by  BF.-\M  -  B.a.se  Layer  I'ixed-Free 
Constrainind  Laver  Free-Free 
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;urr  '21.  Brain  Los.s  Factor  vs.  Thickness  a.s  Predicted  by  BBAM  -  Comparison  of  Cantiii. 
[boundary  Conditions,  Mode  d 
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VIII.  Recommendations  for  Further  Research 


There  are  a  number  of  things  which  should  be  done  to  follow  up  on  these  results.  The  formu¬ 
lation  should  be  generalized  to  more  than  three  layers,  and  the  software  modified  accordingly.  A 
similar  approach  should  be  applied  to  more  complicated  geometries,  and  results  compared  with  pre¬ 
vious  findings  (18:107).  A  parametric  study  should  be  done  to  gain  a  more  thorough  understanding 
of  optimum  design  configurations. 


IX.  Conclusion 


The  equations  of  motion  of  a  three  layer  beam  were  derived  using  Hamilton’s  principle. 
Lateral  displacement  of  the  composite  beam,  and  longitudinal  displacements  of  the  eleistic  layers, 
were  considered  separately.  The  formulation  resulted  in  a  set  of  eight  boundary  conditions  which 
had  to  be  satisfied  in  order  to  calculate  complex  frequencies  for  damped  normal  modes.  The 
boundary  conditions  had  nontrivial  solutions  only  when  the  determinant  of  the  boundary  condition 
matrix  was  zero. 

Tlie  boundary  conditions  matrices  were  presented  for  four  important  cases,  selected  to  il¬ 
lustrate  the  effects  of  permitting  mixed  boundary  conditions  on  frequency  and  loss  factor  for  the 
damped  normal  modes. 

Numerical  results  were  presented  for  a  number  of  test  cases  for  each  of  the  four  boundary 
conditions.  First,  cases  were  run  for  which  simple  theoretical  solutions  were  available,  to  determine 
whether  the  formulation  as  presented  would  correctly  predict  the  frequencies,  damping,  and  mode 
shapes.  The  results  indicated  that  the  software  was  performing  correctly.  Second,  some  cases  were 
shown  to  illustrate  the  effect  of  permitting  the  constraining  layer  to  have  a  different  boundary 
condition  than  the  base  layer.  For  the  cantilever  case,  for  all  but  the  first  mode,  fixing  the  top  layer 
at  X  =  0  resulted  in  higher  damping  than  permitting  it  to  be  free  at  x  =  0.  For  the  simply  supported 
cases,  fixing  the  top  layer  such  that  it  rotates  with  the  base  layer  at  x  =  0  resulted  in  lower  damping 
than  permitting  the  top  layer  to  be  free  at  x  =  0.  Finally,  the  effects  of  considering  mixed  boundary 
conditions  on  optimum  design  was  illustrated  by  calculating  loss  factor  as  a  function  of  adhesive 
layer  thickness  for  all  boundary  conditions.  The  thickne.ss  at  which  peak  damping  occured  shifted 
dramatically  for  some  modes,  and  very  little  or  not  at  all  for  others.  The  results  as  a  whole  suggest 
that  the  optimum  design  changes  depending  on  whether  the  top  layer  has  the  same  boundary 


conditions  as  the  base  layer. 


Appendix  A.  Notation 


Variables  are  shown  in  roughly  the  order  in  which  they  occur  in  the  text.  For  all  quantities, 
the  substript  R  denotes  the  real  part,  and  the  substript  1  denotes  the  imaginary  part.  A  dot 
denotes  differentiation  with  respect  to  time. 

<7  =  stress 

f  =  strain 

E  =  Young’s  Modulus,  experimentally  determined  proportionality  constant  in  the  rela¬ 
tionship  of  a  general,  one  dimensional,  linear  elastic  material 

El  =  An  experimantally  determined  proportionality  constant  the  Kelvin  or  V’oigt  model 
of  viscoelasticity  (3:19) 

£'•_<  =  An  experimantally  determined  proportionality  constant  the  Kelvin  or  Voigt  model 
of  viscoelasticity  (3:19) 

t  =  time 

Qm  and  =  emperical  proportionality  constants  in  a  model  of  a  general,  one  dimensional,  linear 
viscoelastic  material 

M  and  N  =  the  number  of  derivatives  and  proportionality  constants  considered  in  a  particular 
model 

(7q  =  average  stress 

(Q  =  average  strain 

u!  =  complex  frequency 

i  —  square  root  of  -1 

T  =  temperature 

OT  =  reduced  temperature  (17:1-13,503-511) 

A(1)..A(6)=  curve  fit  parameters  in  the  reduced  temperati"'e  equation  (17:1-13,503-  51 1 ),  Eqn  6, 
for  a  given  viscoelastic  material  (17) 

..  _  [j _ L_l^ 
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p  _  _1 _ 1_ 

-  Am  >1(1) 

Cc  =  A{b)  -  ^(4) 

Dc  =  A(6)  -  >1(4) 

De  =  DbCa  —  CbDa 

a  —  {DbCc  —  Cb  Dc)  / De 

b=  {CaDc~  DaCc)IDe 
G  =  Complex  shear  modulus 
/  =  frequency 

i?(l) . . .  D{6)  =  curve  fit  parameters  in  the  shear  modulus  equation  (17:1-13,503-511),  Eqn  7,  for 
a  given  viscoelastic  material  (17) 

X  =  independant  variable  measuring  distance  along  the  length  of  the  beam 
T  —  total  kinetic  energy  of  the  beam 
M  =  total  mass  in  a  general  second  order  system 
u  =  displacement  for  a  general  second  order  system 
k  =  spring  constant  or  eleistic  modulus  in  a  general  second  order  system 
L  =  length  of  the  beam 
6  =  width  of  the  beam 
T)  =  loss  factor 

T)„  =  loss  factor  of  a  viscoelastic  material 
T)(,  =  equivalent  loss  factor  of  the  layered  beam 
Peo  =  density  of  base  (elastic)  layer 
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teo  =  thickness  of  the  base  layer 


t^fO  =  lateral  displacement  of  the  base  layer 

UfO  =  longitudinal  displacement  of  the  base  layer 

pt,i  =  density  of  1st  elaistic  layer 

<1,1  =  thickness  of  the  1st  viscoelastic  layer 

Vt,i  =  lateral  displacement  of  the  1st  viscoelastic  layer 

Ui  i  =  longitudinal  displacement  of  the  1st  viscoelastic  layer 

Pel  =  density  of  1st  elastic  layer 

<ei  =  thickness  of  the  1st  cleistic  layer 

I’ei  =  lateral  displacement  of  the  1st  elastic  layer 

Uei  =  longitudinal  displacement  of  the  1st  elastic  layer 

V  =  lateral  displacement  of  the  composite  beam,  assuming  the  lateral  displacement  of  all 

layers  is  the  same 

(p<)r  =  equivalent  areal  density  (mass  per  unit  area)  of  the  beam 

V  =  total  potential  energy  of  the  beam 

Eeo  =  modulus  of  elasticity  of  the  base  (elastic)  layer 
Gtii  =  shear  modulus  of  the  1st  viscoelastic  layer 
Efi  =  modulus  of  elasticity  of  the  1st  elastic  layer 

di  =  \  +  (<«i  +  <eo)/(2<i.l) 

Dt  =  (E,ot%)/[\2{l  - +  (E,,<3,)/[i2n  _  ul)] 

7,,  I  =  shear  strain  in  the  first  viscoelastic  layer 
=  (u,i  -  U(o)/<t,i  +  d\dv/dx 
S  =  variational  operator 

j  =  index  identifying  the  eight  roots  of  the  characteristic  equation,  Eqn  29,  and  the  con¬ 
stants  associated  with  each  root 
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Pj  —  A  root  of  the  characteristic  equation 
.4j  =  average  displacement  v  associated  with  pj 
Bj  =  average  displacement  Ujo  associated  with  Pj 
Cj  —  average  displacement  Uei  associated  with  pj 
=  B,/A, 

7j  =  CjMj 

=  /?;  -  7; 

1^3  =  -  /^j)  +  d-\UiP}  -  (Dtp])IG^i 

Tj  =  ~ij  +  t^.idiPj 

BEAM  =  name  of  the  computer  program  written  as  a  part  of  this  research  effort,  to  imple¬ 
ment  the  equations  presented  in  the  thesis 

PRED’1’=  name  of  the  program  published  in  the  damping  design  guide  in  Reference  (16:5-29) 
A  =  p-m/(£/)]'/‘'  for  a  homogeneous  cantilever  beam  (10:224-227) 

I  =  moment  of  inertia  of  homogeneous  cantilever  beam 
m  =  mass  per  unit  length  of  homogeneous  cantilever  beam 
n  =  the  mode  number 
D  =  a  constant 
TT  =  3.14159  . . , 

F  =  a  constant 
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Appendix  B.  Description  of  User  Input  for  BEAM 


To  use  the  software  written  as  part  of  this  research  (BEAM),  an  input  data  file  must  be 
prepared.  Following  is  a  description  of  the  format  of  the  input  data  file,  and  an  explanation  of 
input  variables. 

The  first  nine  characters  of  each  input  data  line  are  ignored  by  the  software.  This  makes  it 
possible  to  put  the  name  of  the  variable  on  the  line,  followed  by  the  data,  making  it  easier  to  keep 
track  of  things.  The  data  is  then  read  in  free  format  fror;.  columns  10-80.  The  location  of  the  data 
on  the  line  is  not  important,  but  data  must  be  included  for  every  variable.  If  a  number  is  zero,  it 
is  not  sufficient  to  leave  that  area  blank.  A  zero  must  be  entered.  This  is  the  price  which  must  be 
paid  for  the  freedom  to  locate  the  data  anywhere  on  the  line.  For  data  lines  on  which  more  than 
one  input  variable  must  be  placed,  the  data  must  be  separated  by  normal  field  separators,  such 
as  a  space  or  a  comma.  A  description  of  each  of  the  variables  follows.  For  the  physical  data,  self 
consistent  units  should  be  used. 

OUTPUT  Name  of  the  file  where  output  data  (except  plot  data)  is  to  be  written.  It  may  contain 
a  full  pathname,  i.e.  device:directory/filename 

MESSAG  Name  of  the  file  where  messages  are  to  be  written.  Messages  include  errors,  progress 
reports,  etc.  It  may  contain  a  full  pathname,  i.e.  device;directory/filename 

PLTFIL  Name  of  the  file  where  plot  data  is  to  be  written.  It  may  contain  a  full  pathname,  i.e. 
deviceidirectory/filename 

ICASE  A  number  which  identifies  the  case  being  run.  This  number  is  not  used  by  the  program, 
but  it  is  included  in  the  database  output  described  below. 
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CDESC  A  character  string  describing  the  case  being  run.  It  is  not  used  by  the  program,  but  it 
is  included  in  the  databcise  output  described  below. 

IBC  A  number  indicating  the  boundary  condition  to  be  considered.  Boundary  conditions 
currently  supported  are 

1.  All  layers  fixed  at  i  =  0  and  free  at  i  =  £,. 

2.  All  layers  fixed  at  both  ends. 

3.  Base  layer  fixed  at  a;  =  0  and  free  at  i  =  L.  Constraining  layer  free  at  both  ends. 

4.  Base  layer  pinned  at  a;  =  0  and  a.t  x  =  L.  Constraining  layer  free  at  both  ends. 

5.  Base  layer  pinned  at  i  =  0  and  at  i  =  L.  Constraining  layer  rotates  with  base 
layer  at  x  =  0,  and  is  free  ad  x  =  L. 

BCDESC  A  character  string  describing  the  boundary  condition.  It  is  not  used  by  the  program, 
but  it  is  included  in  the  database  output  described  below. 

PROMPT  A  logical  variable  (T  or  F).  If  true,  the  software  will  ask  prompt  the  user  for  responses, 
such  as  continuing  in  the  event  of  an  error,  or  closing  files,  etc.  If  false,  the  software 
does  not  ask  for  any  user  input.  This  is  useful  for  running  in  batch  mode. 

PLOT  A  logical  variable  (T  or  F).  If  true,  the  software  will  write  the  mode  shapes  as  a  function 
of  position  to  the  file  PLTFIL  in  a  format  which  may  be  read  by  a  plotting  package. 

MODOUT  A  logical  variable  (T  or  F).  If  true,  the  software  will  write  ‘mode’  data,  the  constants 
Ai  . .  .A^,  Bj  and  Bs  to  OUTFIL. 

VERBC  A  logical  variable  (T  or  F).  If  true,  the  software  will  substitute  ‘zeros’  found  back  into 
the  boundary  condition  equations,  calculate  the  error,  and  write  the  error  to  OUTFIL. 

SBOUT  A  logical  variable  (T  or  F).  If  true,  the  software  will  write  result  data  in  a  format  which 
may  be  read  into  a  database  manager.  This  output  only  makes  sense  for  a  databeise 
file  which  is  properly  defined.  The  best  way  to  understand  the  output  is  to  look  at  it 
for  a  couple  of  weeks. 

ZEROUT  A  logical  variable  (T  or  F).  If  true,  the  software  will  write  values  of  complex  frequencies 
at  the  relative  minima. 

UNIFRM  A  logical  variable  (T  or  F).  If  true,  the  grid  used  to  search  for  zeros  will  be  everywhere 
uniformly  spaced.  If  false,  the  grid  size  will  increase  geometrically  from  the  minimum 
to  the  meiximum  values,  resulting  in  a  finer  grid  nearer  the  minimum.  This  should  only 
be  used  if  only  the  first  quadrant  in  the  complex  frequency  plane  is  be.ng  searched. 

WATCH  A  logical  variable  (T  or  F).  If  true,  the  software  will  count  iterations  and  print  progress 
every  IPRINT  iterations.  This  is  not  a  good  idea  if  console  output  is  being  redirected 
to  a  file,  as  in  a  batch  run.  It  is  very  helpful  if  console  is  being  directed  to  the  terminal, 
and  the  run  is  expected  to  take  a  long  time. 
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VARYG  A  logical  variable  (T  or  F).  If  true,  the  complex  shear  modulus  will  be  calculated 
at  each  frequency  value  using  the  temperature  and  emperical  constants  listed  below. 
If  false,  the  temperature  and  emperical  constants  are  all  ignored  except  GVBl  and 
GVB2.  GVBl  is  taken  to  be  the  real  part  of  the  complex  modulus,  and  GVB2  is  taken 
to  be  the  imaginary  part.  They  are  held  constant  for  all  values  of  frequency. 

ITERAT  A  logical  variable  (T  or  F).  If  true,  the  variable  specified  by  ITVAR  will  be  iterated 
according  to  the  values  in  XIT.  Note  that  the  entire  grid  will  be  searched  and  all  zeros 
found  for  every  value  of  ITVAR.  This  could  take  a  while,  so  it  is  a  good  idea  to  know 
how  fast  the  calculations  go  before  attempting  to  iterate  anything. 

LOGINC  A  logical  variable  (T  or  F).  If  false,  the  values  in  XIT  are  used  directly  for  ITVAR. 

If  true,  and  if  ITERAT  is  true,  the  variable  listed  in  ITVAR  will  be  incremented 
‘logarithmically’.  In  other  words,  the  data  points  will  be  equally  spaced  on  a  log  axis. 
The  values  in  XIT  will  be  interpreted  accordingly.  The  first  value  assigned  to  ITVAR 
will  be 

10‘"> 

and  the  maximum  value  assigned  to  ITVAR  will  be 

10^". 

The  increment  value,  Xi  will  be  the  difference  between  the  base  ten  logarithm  of  suc¬ 
cessive  XIT  values.  For  example,  if  the  user  would  like  to  calculate  20  values,  from  10 
to  1000,  equally  spaced  on  a  log  axis,  the  three  values  in  XIT  would  be  1,  .1,  3. 

ITVAR  An  integer  indicating  which  variable  is  to  be  iterated.  The  number  refers  to  the  physical 
data  listed  below.  XL  is  variable  1,  EEOR  is  variable  2,  etc.  down  to  XNUEl  which  is 
variable  26. 

XIT  Three  values  controlling  the  iteration  of  ITVAR.  The  first  number  is  the  minimum 
value  (^i),  the  second  the  increment  (r,),  and  the  third  the  maximum  (im)- 

IPRINT  An  integer  specifying  how  often  progress  is  to  be  printed  to  the  console,  if  WATCH 
is  true.  For  example,  if  IPRINT  =  1000,  the  software  will  print  a  progress  message 
every  1000  iterations.  This  is  an  important  safety  valve.  There  are  times  when  the 
user  will  want  to  interrupt  program  execution.  For  example,  it  is  easy  to  inadvertantly 
ask  the  program  to  search  search  a  grid  which  would  take  it  several  months.  There  are 
also  times  when  it  is  clear  that  things  are  not  running  correctly.  In  both  MS-DOS  and 
AmigaDOS,  a  control-C  interrupt  only  takes  effect  when  the  program  is  performing 
terminal  I/O.  Specifying  a  reasonable  value  for  IPRINT  will  force  the  program  to  ‘come 
up  for  air’  periodically,  and  give  the  user  a  chance  to  interrupt  it  gracefully. 

OMO  Two  real  numbers  specifying  the  starting  value  for  the  complex  frequency  grid.  The 
first  number  is  the  real  value,  the  second  the  imaginary. 

OMMAX  Two  real  numbers  specifying  the  ending  value  for  the  complex  frequency  grid.  The 
first  number  is  the  real  value,  the  second  the  imaginary. 

NOMEGA  Two  integers  specifying  how  many  points  to  calculate  in  the  complex  frequency  grid. 

The  first  number  specifies  how  many  real  values  to  calculate,  the  second  how  many 
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imaginary  values  to  calculate.  In  the  current  implementation,  the  number  of  imaginary 
points  is  limited  to  a  maximum  of  2500,  and  an  unlimited  number  of  real  points  may 
be  calculated.  It  would  be  wise  to  experiment  with  a  few  calculations  before  asking 
for  a  large  grid,  to  get  a  feel  for  how  fast  the  computer  will  work. 

DETZER  A  real  number  controlling  the  end  of  a  search  for  a  zero.  If  the  magnitude  of  the 
determinant  at  a  zero  is  less  than  DETZER,  the  search  terminates. 

OMTOL  A  real  number  controlling  the  end  of  a  search  for  a  zero.  If  the  distance  between  two 
successive  frequency  values  is  less  than  OMTOL,  the  search  terminates. 

NMODES  An  integer  specifying  how  many  zeros  are  to  be  printed  in  the  mode  shape  plot  file, 
PLTFIL.  This  does  not  limit  the  number  of  zeros  which  will  be  found,  only  the  number 
which  will  be  written  to  the  plot  file. 

NX  For  the  mode  shape  plot  file,  this  is  the  number  of  X  values  for  which  the  displacements 
will  b''  calculated  and  written.  They  will  be  uniformly  spaced  from  x  =  0  to  x  =  L. 

XL  The  length  of  the  beam. 

EEOR  The  real  part  of  the  modulus  of  elasticity  of  the  base  layer. 

EEOl  The  imaginary  part  of  the  modulus  of  elasticity  of  the  base  layer. 

TEO  The  thickness  of  the  base  layer. 

RHOEO  The  density  of  the  base  layer. 

XNUEO  Poisson’s  ratio  for  the  base  layer. 

GVAl  Emperical  constant  A1  used  to  calculate  the  shear  modulus. 

GVA2  Emperical  constant  A2  used  to  calculate  the  shear  modulus. 

GVA3  Emperical  constant  A3  used  to  calculate  the  shear  modulus. 

GVA4  Emperical  constant  A4  used  to  calculate  the  shear  modulus. 

GVA5  Emperical  constant  A5  used  to  calculate  the  shear  modulus. 

GVA6  Emperical  constant  A6  used  to  calculate  the  shear  modulus. 

GVBl  Emperical  constant  B1  used  to  calculate  the  shear  modulus. 

GVB2  Emperical  constant  B2  used  to  calculate  the  shear  modulus, 

GVB3  Emperical  constant  B3  used  to  calculate  the  shear  modulus. 
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GVB4  Emperical  constant  B4  used  to  calculate  the  shear  modulus. 

GVB5  Emperical  constant  B5  used  to  calculate  the  shear  modulus, 

GVB6  Emperical  constant  B6  used  to  calculate  the  shear  modulus. 

TEMPER  The  temperature,  used  to  calculate  the  shear  modulus. 

TVl  Thickness  of  the  adhesive  layer. 

RHOVl  Density  of  the  adhesive  layer. 

EEIR  The  real  part  of  the  modulus  of  elasticity  of  the  constraining  layer. 

EElI  The  imaginary  part  of  the  modulus  of  elasticity  of  the  constraining  layer. 
TEl  The  thickness  of  the  constraining  layer. 

RHOEl  The  density  of  the  constraining  layer. 

XNUEl  Poisson’s  ratio  for  the  constraining  layer. 

A  sample  input  data  file  follows.  Note  that  the  actual  data  begins  in  column  ten. 


OUTPUT  =  data; test . out 
MESSAG  =  ♦ 

PLTFIL  =  datartest .pit 
ICASE  =  3 

CDESC  =  dcunped-beeuo 

IBC  =  1 

BCDESC  =  fixed/free 

PROMPT  =  F 

PLOT  =  F 

MODOUT  =  T 

VERBC  =  T 

SBOUT  =  F 

ZEROUT  =  F 

UHIFRM  =  T 

WATCH  =  T 

VARYG  =  T 

ITERAT  =  F 

LOGIIC  =  F 

ITVAR  =  20 

XIT  =  0.001  0.001S4 

IPRIHT  =  1000 

OMO  =  10.  -10. 


0.003 
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OMMAX 

= 

2800. 

HOMEGA 

= 

25 

DETZER 

= 

l.D-20 

OMTOL 

= 

l.D-14 

HMODES 

= 

20 

HX 

= 

100 

XL 

= 

.254 

EEOR 

= 

6.89E10 

EEOI 

= 

0. 

TEO 

= 

.00508 

RHOEO 

= 

2770.0 

XMUEO 

= 

.33 

GVAl 

= 

335. 

GVA2 

= 

280. 

GVA3 

= 

390. 

GVA4 

= 

.7E-1 

GVA5 

= 

.  1142 

GVA6 

= 

.3E-1 

GVBl 

= 

.2E6 

GVB2 

= 

1200. E6 

GVB3 

= 

1 .885E7 

GVB4 

.55 

GVB5 

= 

1.5 

GVB6 

= 

.  1 

TEMPER 

= 

350. 

TVl 

= 

.000254 

RHOVl 

= 

969.0 

EEIR 

= 

6.89E10 

EEII 

0. 

TEl 

= 

.000254 

RHOEl 

2770.0 

XNUEl 

= 

.33 

10. 

5 
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Appendix  C.  Programming  Notes  for  Program  BEAM 


The  software  was  written  in  ANSI  standard  Fortran  77,  and  corn^  'led  and  run  on  a  Com¬ 
modore  Amiga  2500/020  using  ABSOFT  Fortran  (1).  The  code  should  be  easily  portable  to  other 
computers  with  only  minor  modifications.  Nonstandard  code  generally  involved  such  things  cis 
system  time  calls,  and  device  or  filename  conventions.  Anywhere  system  dependant  code  was  used, 
it  WcLS  noted  in  the  source. 

Most  of  the  variables  are  declared  in  separate  files  and  placed  in  common  blocks.  The  common 
files  are  then  included  in  each  routine  which  needs  access  to  those  variables.  This  simplifies  making 
changes  to  common  blocks.  Syntax  of  the  INCLUDE  statment  may  vary  among  compilers,  as  will 
the  filename  conventions.  The  include  files  also  have  descriptions  of  the  variables.  All  variables  in 
COMMO.N  are  initialized  in  a  BLOCK  DAT.A  subroutine,  at  the  end  of  the  main  source  file. 

Jo  maximize  portability,  most  variables  were  kept  under  six  characters.  Case  is  significant  in 
ih.'  variables  in  the  system  time  routine,  which  is  machine  specific.  Accordingly,  the  compiler  was 
require  '  to  consider  case  significance.  I'o  avoid  any  problems,  all  code  was  written  in  upper  case. 

The  data  was  collected  and  maintained  in  a  database.  Ad  input  and  all  output  data  for 
each  ‘zero’  located  in  aii  test  cases  was  kept.  The  database  manager  used  ■■  -  Superbase  Profes¬ 
sional  3.0  (12).  but  the  data  can  be  made  available  in  any  of  a  number  of  common  formats.  All 
the  data,  uncomi  ressed,  occupies  approximately  C  megabytes  of  data.  When  exported  to  dBa- 
selll  for.nat  and  compressed  using  LHarc,  the  data  occupies  two  compressed  files  of  approximately 
500  kilobytes  each.  LHarc  is  a  standard,  freely  redistributable  compression  routine  widely  available 
for  most  persona!  computers  and  mainframes. 
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